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1.0  SUMMARY 


In  response  to  a  request  from  Dr.  Jim  Scofield  of  RZPE,  the  RZPS  thermal  group  constructed  an 
apparatus  to  simulate  the  cooling  of  a  silicon  carbide  (SiC)  MOSFET  using  a  mini-channel  heat 
exchanger  with  a  single  phase  (liquid)  working  fluid.  The  chip  was  simulated  with  a  sample  of 
SiC  that  is  heated  on  one  surface  and  cooled  on  the  opposing  surface  using  a  mini-channel  heat 
exchanger.  One  goal  of  this  work  was  to  demonstrate  that  the  apparatus  has  the  capability  to 
remove  a  heat  flux  of  lOOW/cm^  from  the  surrogate  chip  while  maintaining  its  temperature 
below  lOO^C.  A  second  goal  was  to  characterize  the  thermal  and  flow  behavior  of  the  mini¬ 
channel  heat  exchanger  for  different  flow  rates  and  heat  fluxes  by  experiments  and  numerical 
simulation  using  computational  fluid  dynamics.  This  work  is  part  of  a  longer  term  Air  Force 
program  which  has  multiple  projects  seeking  to  find  the  optimum  electronic  module  package. 

The  fabricated  test  stand  and  heat  exchanger  showed  relatively  good  effectiveness  in  transferring 
thermal  energy.  In  the  experiments,  the  maximum  heat  flux  removed  by  the  mini-channel  heat 
exchanger  was  110-117  W/cm^  with  a  maximum  stem  temperature  of  1 12°C  at  the  interface.  The 
Nud  values  determined  here  were  in  the  range  of  9-23,  which  followed  several  heat  transfer 
correlations  reasonably  well.  For  flows  up  to  a  Rcd  of  3600,  there  was  good  agreement  with  the 
theoretical  pressure  loss.  The  transition  to  turbulent  flow  at  a  Rcd  of  3600  for  the  present 
experiments  is  greater  (by  1000)  than  many  other  works  reviewed  and  is  likely  due  to  the  very 
low  surface  roughness  of  the  channels  (0.2  qm).  Easily,  computational  fluid  dynamics 
simulations  provided  insight  into  the  effect  of  the  heated  length  on  Nusselt  number  estimates  and 
flow  and  convective  heat  transfer  behavior. 
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2.0  INTRODUCTION 


The  demand  for  higher  flight  speeds,  more  power-dense  avionics,  and  advanced  weaponry  is 
limited  by  aircraft  weight  and  thermal  restrictions.  The  Air  Force  More  Electric  Aircraft 
Initiative  has  the  goal  of  using  improved  electronic  technology  to  replace  heavy  mechanical 
systems  which  tax  the  energy  generation  of  power  plants  and  may  have  relatively  high 
maintenance  costs.  For  example,  using  electrical  actuation  rather  than  hydraulic  actuation 
eliminates  complex  piping,  mechanical  components,  and  the  need  for  hydraulic  fluid.  However, 
the  replacement  of  traditional  aircraft  mechanical  systems  which  use  well-established  cooling 
methods  with  electrical  systems  which  use  advanced  cooling  methods  and  materials  are  active 
areas  of  research.  One  approach  to  the  thermal  management  of  electronic  circuits  is  to  increase 
the  allowable  operating  temperature  by  the  development  of  new  materials,  such  as  SiC. 


The  use  of  SiC  to  fabricate  electronic  components  offers  advantages  over  the  use  of  Si  which  has 
been  used  in  electronic  components  for  decades.  Moreover,  new  aircraft  electronic  components 
will  have  more  stringent  performance  requirements  than  those  of  the  past.  SiC  offers  higher 
temperature  limits  (200-250°C)  for  electronic  components  relative  to  those  which  use  Si 
(temperature  limits  in  the  range  125-150‘’C).  In  addition,  the  upper  temperature  limit  for  SiC 
components  is  expected  to  reach  600“C  (for  wind  turbines  and  hybrid  energy  generation)  which 
will  drive  future  research  and  development  (Marckx,  2006).  In  the  use  of  high-voltage  electronic 
devices,  the  breakdown  voltage  (voltage  at  which  electrons  begin  unwanted  travel  across 
interfaces)  is  an  important  characteristic  to  consider.  SiC  has  a  much  higher  breakdown  voltage 
(-'3MV/cm)  than  Si  (0.6  MV/cm)  which  allows  SiC  devices  to  have  relatively  high  operating 
voltages  and  small  dimensions  without  power  leakage  to  surrounding  elements  (Ozpineci,  2007). 
Small  devices  are  often  necessary  for  space-constrained  environments.  Currently,  only  a  few  SiC 
devices  exist  which  include  the  junction  field  effect  transistor  (JFET)  and  Schottky  diode. 
However,  these  are  hybrid  devices  comprised  of  Si  insulated  gate  bipolar  transistors  (IGBT"s) 
and  SiC  switching  components.  Metal  oxide  semiconductor  field  excited  transistor  (MOSFET) 
devices  are  the  future  of  SiC  devices  with  release  possible  for  the  near  future  (Hasanuzzaman, 
2004).  MOSFET  devices  do  not  require  a  voltage  to  switch  off  (unlike  a  JFET),  and  high 
frequency  capabilities  (greater  than  100  kHz)  are  capable  with  low  internal  capacitance  which 
reduces  full  on-off  time  and  electrical  loss  (Bontemps,  2007). 


Unfortunately,  the  use  of  SiC  as  a  MOSFET  device  has  drawbacks  because  of  technology  limits 
with  cooling  and  power  connections.  Presently,  the  use  of  SiC  in  electronic  devices  is  limited  to 
200°C  because  of  the  oxidation  problem  with  the  electrical  contacts.  The  oxidation  problem 
arises  from  the  requirements  of  the  contacts  to  be  Ohmic  (they  contain  internal  resistance)  and 
the  method  of  gate  attachment  (Carrion,  2007).  At  temperatures  above  200°C,  the  dielectric  layer 
(SiOi)  which  attaches  the  gate  to  the  device  becomes  unstable  and  increases  the  material 
resistance. 


The  relatively  large  voltages  and  high  operation  frequencies  (on  and  off  switching)  that  are 
expected  with  the  use  of  SiC  devices  translate  to  large  levels  of  waste  heat  production  per  area. 
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The  continuing  trend  of  decreasing  external  surface  area  results  in  a  high  heat  flux  level  which 
exacerbates  thermal  management  problems  and  leads  to  higher  electronic  device  temperatures. 


Keeping  the  operation  temperature  of  SiC  devices  below  200°C  in  order  to  take  advantage  of 
high  voltage  and  high  speed  capabilities  is  a  significant  task.  To  meet  this  task,  work  is  needed  in 
the  design  and  development  of  material  properties  and  cooling  systems.  Standard  modules  use 
DBC  (direct  bonded  copper,  copper-AlN-copper)  as  the  intermediate  between  the  heat  sink  and 
electronic  device.  However  with  the  higher  temperature  of  operation  expected  with  SiC, 
verification  must  be  made  that  DBC  can  still  perform  properly.  The  combined  properties  of 
coefficient  of  thermal  expansion  (CTE),  thermal  conductivity,  and  dielectric  value  for  the 
substrate  materials,  solder,  and  the  electronic  device  will  have  to  meet  more  closely  the  demands 
of  evolving  SiC  devices  (Griffin  and  Koebe,  2006).  Griffin  and  Koebe  (2006)  performed  testing 
where  a  module  of  SiC  with  copper-AlN-copper  successfully  operated  up  to  200°C,  but  failures 
plagued  their  testing  when  temperature  cycling  became  too  rapid.  However,  their  work  still 
shows  that  there  is  future  potential  with  SiC  devices  and  that  more  research  is  needed  in  this 
area. 


Encouraging  future  testing  is  being  planned  by  the  NREE  for  wind  turbine  applications. 
Researchers  there  have  performed  modeling  experiments  for  three  configurations  of  materials, 
the  most  notable  set  being  AlSiC  and  AIN.  A  SiC  device  operating  with  4160  VAC,  208  Amps, 
at  3  kHz  was  shown  to  operate  without  thermal  stress  failures  at  300“C.  However,  fluctuating 
temperature  (cycling)  testing  was  not  reported  (Marckx,  2006).  The  fabrication  of  electronic 
devices  that  will  not  fail  by  material  and  joint  fatigue  under  high  cycling  rates  and  temperatures 
is  a  significant  challenge.  This  may  require  material  and  solder  changes  or  more  advanced 
layering  schemes.  As  temperatures  and  cycling  rates  increase,  the  mismatch  of  CTE  and 
expansion  rates  for  each  module  material  becomes  critical  as  the  de-bonding  of  materials 
becomes  more  likely  (Marckx,  2006). 


Many  research  groups  have  performed  tests  (Griffin  and  Koebe,  2006;  Mossor,  1999;  Schulz- 
Harder  et  ah,  2007;  and  Hasanuzzaman  et  ah,  2004)  with  various  materials  for  high  temperature 
operation,  but  comparison  of  the  results  is  difficult  due  to  differences  in  testing  setup,  type  of 
electronic  device,  and  range  of  operating  parameters.  Thus,  there  is  a  need  for  more  research 
involving  the  experimental  and  numerical  simulation  of  electronic  devices,  especially  with 
regard  to  material  behavior  and  thermal  management.  Performing  research  on  electronic  device 
materials  presents  an  opportunity  for  testing  with  mini-  or  micro-  channels  as  the  heat  sink.  Mini- 
and  micro-  channels  are  reported  to  have  significant  capabilities  for  removing  high  heat  fluxes. 


In  response  to  a  request  from  Dr.  Jim  Scofield  of  RZPE,  the  RZPS  thermal  group  has  constructed 
an  apparatus  to  simulate  the  cooling  of  a  SiC  MOSEET  using  a  mini-channel  heat  exchanger 
with  a  single  phase  (liquid)  working  fluid.  The  chip  was  simulated  with  a  piece  of  SiC  that  is 
heated  on  one  surface  and  cooled  on  the  other  using  a  mini-channel  heat  exchanger.  One  goal  of 
this  work  was  to  demonstrate  that  the  apparatus  has  the  capability  to  remove  a  heat  flux  of 
lOOW/cm^  from  the  surrogate  chip  while  maintaining  its  temperature  below  200°C.  A  second 
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goal  was  to  characterize  the  thermal  and  flow  behavior  of  the  mini-ehannel  heat  exehanger  for 
different  flow  rates  and  heat  fluxes  by  experiments  and  numerical  simulation  using 
computational  fluid  dynamics.  This  work  is  part  of  a  longer  term  Air  Force  program  whieh  has 
multiple  projects  seeking  to  find  the  optimum  eleetronic  module  paekage. 
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3.0  EXPERIMENTAL  DESCRIPTION 


A  test  apparatus  was  designed  to  provide  a  way  to  determine  the  heat  transfer  through  a  SiC 
sample  and  supporting  materials.  (In  the  future,  other  materials  and  material  stacks  may  be  used.) 
These  materials  were  cooled  with  a  mini-channel  heat  exchanger  that  was  designed  and  tested  as 
a  part  of  this  project.  Figure  1  is  a  block  diagram  that  shows  the  experimental  components  of  the 
test  cell  (test  stand  and  required  components).  Figure  1  indicates  that  a  cooling  fluid  flows  from 
the  cooler/pump  and  through  a  40  micron  filter  before  reaching  the  mini-channel  heat  exchanger 
(within  the  test  stand  enclosure).  Distilled  water  was  used  to  avoid  contaminant  materials  that 
could  damage  the  flow  meter  or  interior  passages  of  the  chiller.  The  heat  exchanger  transfers 
thermal  energy  from  a  source  (copper  heater  and  stem)  to  the  fluid.  After  flowing  through  the 
heat  exchanger,  the  water  passes  through  the  flow  meter.  The  flow  meter  contains  a  magnetic 
turbine  and  is  calibrated  to  provide  the  mass  flow  rate  (Appendix  3.3).  The  last  components 
before  the  water  returns  to  the  cooler  are  the  two  flow  valves.  The  valves  are  used  in  parallel  for 
small  flow  adjustments.  Each  of  the  components  and  supporting  instrumentation  of  Figure  1  are 
described  in  detail  in  following  paragraphs.  In  addition,  an  image  of  the  test  components  and 
supporting  instrumentation  (data  acquisition  (DAQ)  and  power  supplies)  are  shown  in  Figure  2 
on  the  supporting  cart  fabricated  for  transport  of  the  apparatus. 
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Figure  1.  Schematic  of  Test  Cell 
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Figure  2.  Experimental  Electronics  Cooling  Apparatus 


3.1  Test  Stand 

The  test  stand  was  the  primary  group  of  experimental  components  for  this  project.  „Test  stand" 
refers  to  the  mini-channel  heat  exchanger,  the  copper  block  (contains  an  electrical  heater),  and 
the  aligning  fixtures.  Two  test  stand  designs,  the  second  was  an  evolution  of  the  first,  were 
completed  for  the  testing.  The  first  design  (Figure  3)  was  smaller  and  involved  fewer 
components  than  the  second.  The  heat  exchanger  in  the  first  design  had  channel  lengths  of  '-1 1 
mm.  Unfortunately,  this  initial  design  had  repeatability  problems  with  the  experiments  and 
required  re-design  of  the  test  stand.  The  heat  exchanger  was  awkwardly  bent  which  made  fiat 
alignment  with  the  heating  surface  impossible.  Also,  the  corner  supports  allowed  the  heat 
exchanger  to  slide. 
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Figure  3.  The  First  Test  Stand 


The  design  of  the  seeond  test  stand  had  three  goals  set  to  insure  a  better  outeome.  First,  the  heat 
exchanger  needed  to  be  firmly  held  in  place  above  the  copper  stem.  The  heater  (electric)  was 
contained  within  the  copper  heater  block  shown  in  Figure  3.  At  the  top  of  the  heater  block,  the 
material  converges  to  a  smaller  square  shape,  the  stem.  The  stem  either  by  direct  contact  with  the 
heat  exchanger  or  through  a  material  sample  transfers  heat  to  the  fluid.  Thus,  good  contact  and 
positioning  are  important.  The  second  goal  was  to  have  the  heat  exchanger  alignment  with  the 
heater  be  adjustable.  More  specifically,  the  ability  to  adjust  the  position  of  stem  contact  under  the 
heat  exchanger  was  desired.  The  last  goal  was  that  the  test  stand  needed  to  protect  the  heat 
exchanger  from  external  convective  air  cooling  and  thermal  conduction  to  the  supporting  cart. 


Figure  4  shows  an  image  of  the  second  stand  which  will  be  described  in  the  remainder  of  the 
experimental  description.  Figure  4  shows  the  copper  heater,  convection  shield,  the  alignment 
structures,  and  a  pressure  transducer.  The  heat  exchanger  is  located  at  the  center  of  the  test  stand. 
Figure  4  shows  a  rigid  alignment  stand  that  was  made  for  the  second  design.  The  second  test 
stand  provided  better  repeatability  in  the  experimental  measurements  due  to  more  rigid  and 
better  fabricated  positioning  parts.  A  detailed  AutoCAD  drawing  is  shown  in  Figure  5. 
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Figure  4.  The  Second  Test  Stand 
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Figure  5.  Section  Drawing  of  the  Test  Stand 


3.1.1  Copper  Heater  Bloek 

The  copper  heater  block  was  designed  to  contain  a  cylindrical  resistance  heater  (3/8  in  diameter 
X  2  in  length)  and  converge  to  a  small  (8  mm  x  8  mm)  square  where  the  samples  would  be 
located.  To  promote  a  one-dimensional  temperature  distribution  in  the  stem,  high  thermal 
conductivity  copper  was  selected  as  the  material.  Figure  6  shows  details  of  the  copper  block 
which  holds  the  electrical  heater  (3/8  in  diameter  and  --2  in  long).  The  stem  which  comes  in 
contact  with  the  heat  exchanger  is  located  at  the  top  of  the  copper  block. 
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Figure  6.  Sectioned  Copper  Heater  Block 


3.1.2  The  Heat  Exchanger 

Heat  sinks  and  heat  exchangers  are  commonly  used  for  high  temperature  electronic  applications. 
As  operating  conditions  become  more  difficult,  proper  design  will  become  more  critical.  Mini¬ 
channel  heat  exchangers  (rectangular  or  round)  have  proven  to  be  efficient  at  transferring  large 
heat  fluxes.  The  channels  for  the  present  mini-channel  heat  exchanger  (635  pm  round  channels 
in  rectangular  nickel  substrate)  were  manufactured  under  the  direction  of  Grady  Yoder  at  Oak 
Ridge  National  Laboratory  by  electro-plating  nickel  on  parallel  fibers.  The  fibers  were  then 
removed  by  acid,  resulting  in  channels  that  were  slightly  unparallel  and  out  of  center,  but 
smooth.  The  expected  surface  roughness  for  this  fabrication  technique  is  15nm  (Papautsky  et  ah, 
1999). 


Figure  7a  shows  an  image  of  the  mini-channel  heat  exchanger  with  the  inlet/exit  stainless-steel 
tubes.  To  simplify  the  fluid  dynamics,  it  is  desirable  to  have  the  tubing  that  connects  with  the 
heat  exchanger  to  be  in  a  straight  line  with  the  heat  exchanger  channels  with  few  bends. 
However,  the  design  also  needs  to  accommodate  secure  mounting  over  the  stem  and  provide 
good  alignment.  An  image  of  the  mini-channel  heat  exchanger  is  shown  in  Figure  7a.  Figure  7b 
shows  the  three  parts  of  the  mini-channel  heat  exchanger:  the  mini-channel  section  (green), 
mixing  chambers  (blue,  purple),  and  the  inlet  and  exit  tubes  (yellow). 
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(a) 


(b) 

Figure  7.  An  Image  (a)  and  a  Drawing  (b)  of  the  Mini- Channel  Heat  Exchanger 


2 

The  heat  exehanger,  required  to  transfer  100  W/em  ,  needed  to  transfer  heat  through  a  small  area 
(8  mm  X  8  mm).  For  the  eharaeterization  of  the  flow  and  heat  transfer  within  the  ehannels,  it  was 
not  necessary  to  use  material  samples  for  all  tests.  That  is,  the  copper  stem  rather  than  a  SiC 
sample  was  in  contact  with  the  heat  exchanger  for  most  of  the  testing.  The  contacting  area  of  the 
stem  and  the  heat  exchanger  were  the  same  (8  mm  x  8  mm)  such  that  the  heat  exchanger  mates 
with  the  stem.  One  of  the  methods  to  determine  the  heat  input  to  the  heat  exchanger  was 
measurement  of  the  temperature  increase  of  the  cooling  water.  The  heat  exchanger  width  was 
chosen  to  closely  match  the  stem  size.  The  total  heat  exchanger  width  of  10.54  mm  corresponds 
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to  1 1  circular  channels  with  a  diameter  of  635  |um  and  spaced  323  pm  eenter-to-center  (Figure 
8).  One  less  channel  could  have  been  used,  but  room  was  left  for  small  misalignment. 


The  heat  exchanger  ehannels  are  single  pass  and  parallel  to  eaeh  other.  The  delivered  thiekness 
was  1.778  mm  with  the  channels  located  offset  from  the  center  of  the  material,  nickel.  The  heat 
exchanger  channels  were  used  with  the  thicker  side  down  to  allow  for  machining.  Figure  8a 
(dimensions  are  in  mm)  shows  the  heat  exehanger  channels  without  the  attached  accessories  of 
Figures  7a  and  b.  The  eutout,  indicated  as  8  mm,  is  aetually  wider  to  allow  for  easier  positioning 
of  the  stem  with  the  heat  exchanger.  All  of  the  bottom  material  could  have  been  removed  but  was 
left  for  heat  spreading  affects. 


(a) 


(b) 

Figure  8.  Angled  (a)  and  Front  (b)  View  of  the  Mini- Channel  Heat  Exchanger 


Design  of  the  heat  exchanger  also  involved  determining  the  optimum  channel  length.  Although 
turbulent  flow  generally  results  in  higher  heat  transfer  rates  than  with  laminar  flow,  the  pressure 
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loss  with  laminar  flow  in  this  channel  is  lower  than  that  for  turbulent  flow.  Thus,  having  laminar 
flow  with  developing  hydrodynamie  boundary  and  thermal  layers  along  the  length  of  the 
ehannels  was  a  goal.  If  the  ehannel  length  was  too  long  and  both  the  temperature  and  veloeity 
profiles  were  fully  developed,  the  eooling  would  not  be  optimal.  In  addition,  the  laminar  flow 
may  transition  to  turbulent  flow.  The  ehannels  were  designed  such  that  the  heated  section  was 
located  before  fully  developed  hydrodynamie  eonditions  for  most  Reynolds  numbers  (Reo) 
values  eonsidered.  For  the  635|um  diameter  channels  (as  given  by  Eqn.  1),  the  fully  developed 
hydrodynamie  length  is  -“O.Ol  Im  at  a  Reo  of  350  (which  corresponds  to  the  lowest  flow  rate)  and 
is  0.095m  at  a  Reo  of  3000.  (Wileox,  2000) 

|  =  0.05*Re^.  (1) 

The  actual  heat  exehanger  ehannel  lengths  were  21.59  mm  which  locates  the  region  of  a  fully 
developed  veloeity  profile  at  the  mid-point  of  the  channels. 


In  fabrieating  the  heat  exchanger,  a  method  of  attaehing  the  inlet  and  exit  flow  tubes  (stainless- 
steel  625)  to  the  ehannels  had  to  be  considered.  In  addition,  the  flow  must  be  well-mixed  after 
leaving  the  ehannels  so  that  thermoeouple  measurements  eould  reasonably  represent  the  bulk 
water  temperature  inerease.  A  mixing  ehamber  design  (Figure  9)  was  seleeted  for  attaehment  of 
the  heat  exehanger  ehannels  and  the  inlet/exit  tubing.  The  mixing  ehambers  (Ineonel  600)  were 
maehined  as  two  pieees  whieh  were  later  welded  together  to  make  one  ehamber.  The  purple  half 
of  the  mixing  ehamber  (Figure  9)  has  a  reetangular  slot  where  the  heat  exchanger  was  inserted 
and  attaehed  (Laser-welded  -  Preeision  Joining  Teehnologies).  The  other  half  eontains  a  hole  for 
the  inlet  or  exit  fluid  tube.  There  is  an  8:1  ratio  in  eross-seetional  area  ehange  between  the 
mixing  ehambers  and  the  1 1  channels  to  minimize  pressure  loss  and  to  provide  even  flow  into 
the  channels.  The  mixing  chamber  internal  dimensions  are  1.5  mm  x  8.382  mm  x  12.2  mm.  The 
stainless-steel  inlet/exit  tubes  (yellow)  have  an  area  ratio  ~  1:1  with  the  mixing  chambers  to 
minimize  the  pressure  loss  at  the  attaehment  region. 
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Figure  9.  Two  Halves  of  the  Mixing  Chamber  and  the  Inlet  (or  Exit)  of  the  Flow  Tube 


3.1.3  Alignment  Stand 

The  alignment  stand  holds  the  heat  exchanger  and  the  copper  stem  in  specific  positions  relative 
to  each  other.  The  alignment  stand  must  be  able  to  allow  the  heat  exchanger  assembly  to  move 
small  distances  for  installation  and  material  expansion  when  heated.  The  stand  also  should  have 
a  way  for  applying  a  load  to  the  top  surface  of  the  heat  exchanger  to  ensure  good  contact 
between  the  heat  exchanger  and  the  stem.  Figure  5  shows  the  assembled  alignment  stand  (the 
upper  plate,  legs,  and  table).  The  mounting  table  contains  the  aligning  holes  (optical  breadboard) 
used  to  locate  the  upper  plate  to  the  heater.  The  upper  plate  (Figure  10)  was  the  alignment 
component  for  the  heat  exchanger;  the  heater  and  stem  were  aligned  with  an  alignment  block 
which  squarely  located  them  with  the  test  stand  legs.  Provisions  for  the  heat  exchanger  and  the 
legs  are  indicated  in  Figure  10b. 
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(a) 


Figure  10.  Top  (a)  and  Bottom  (b)  Views  of  the  Upper  Alignment  Plate 

The  upper  alignment  plate  fabrieated  from  aluminum  had  notehed  loeations  (1  in  x  1  in  x  0.25  in) 
whieh  mated  with  the  supporting  legs  and  a  eenter  seetion  removed  for  the  Glas-Therm/heat 
exehanger  assembly  (Figure  11a).  Positioning  the  heated  stem  to  have  thermal  eontact  with  the 
heat  exehanger  was  important  for  proper  distribution  of  the  applied  heat  flux.  An  alignment 
bloek  (Figure  1  lb)  made  of  phenolie  plastie  was  used  to  loeate  the  heater  at  the  eenter  of  the  test 
stand  using  the  referenee  leg.  The  heater  stand  was  made  of  a  few  materials  (Maeor,  aluminum, 
and  phenolie  plastic)  for  insulating  purposes,  vertical  adjustment,  and  stable  positioning. 
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3.1.4  Glas-Therm  Block 


The  heat  exchanger  required  a  block  to  position  it  in  the  upper  alignment  plate  which  would 
allow  for  vertical  movement  and  minimize  heat  loss.  A  heat  exchanger  mounting  block  was 
fabricated  from  Glas-Therm  which  has  thermal  conductivity  value  of  0.6  W/mK  and  was  easy  to 
machine.  The  block  was  rectangular  with  a  heat  exchanger  chamber  cutout  (Figures  11a  and  b). 
The  heat  exchanger  was  held  in  the  block  using  a  ceramic  epoxy  (Omegabond  1500  °F). 


(a) 


(b) 

Figure  11.  Glas-Therm  Block  with  Respect  to  Tuhing  (a)  and  Upper  Alignment  Plate  (b) 
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3 . 1 . 5  Load  Apparatus 

For  enhanced  experimental  repeatability,  a  load  apparatus  was  created  to  apply  and  record  a 
force  on  the  heat  exchanger  and  copper  stem.  This  force  translates  to  a  uniform  contact  pressure 
between  the  heat  exchanger  and  the  stem.  The  load  cell  is  located  at  the  top  of  the  system 
(mounts  on  top  of  upper  plate,  Figures  5  and  12).  Components  of  the  load  apparatus  include  the 
load  screw,  spring  seats,  plunger,  spring,  and  aluminum  block.  The  knob  (Figure  12)  (3/8  in  fine 
thread)  was  used  for  making  the  pressure  adjustment.  The  knob  was  threaded  to  the  upper  spring 
seat  and  applied  a  load  through  a  spring  (spring  constant  of  100  Ibf/in)  to  the  plunger.  The  system 
was  fully  seated  (spring  seat  to  spring  seat)  when  the  system  was  at  steady  conditions.  This 
prevented  overloading  of  the  load  cell  during  warm-up  periods.  The  aluminum  plate  was  for 
mounting  the  knob  and  securing  the  components  to  the  upper  plate.  A  stainless-steel  tube  held 
the  aluminum  plate  and  knob.  The  load  was  applied  to  the  top  surface  of  the  Glas-Therm  block. 
The  load  was  transferred  to  the  heat  exchanger  via  the  load  cell  and  an  aluminum  plunger.  The 
plunger  had  a  1  Lt  in  diameter,  approximately  the  same  width  as  the  Glas-Therm  block. 
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Figure  12.  Apparatus  for  Imposing  Repeatable  Thermal  Load  on  the  Copper  Stem  and 

Heat  Exchanger 
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3.1.6  Test  Stand  Insulation  and  Heat  Loss  Reduction 

To  limit  the  convective  heat  transfer  from  the  copper  block,  a  convection  shield  was  fabricated  to 
surround  the  test  stand.  The  shield  was  mounted  to  the  outside  of  the  test  stand  legs  creating  a 
small  chamber  around  the  heat  exchanger.  Insulating  the  components  was  done  to  limit  transfer 
to  the  heat  exchanger.  Insulating  the  copper  block  (stem  and  heater  portions)  required  many 
attempts  because  the  thermocouples  that  were  attached  to  the  stem  frequently  broke.  Only  one 
layer  (2.54  mm)  of  insulation  (Unifrax,  550K-970J)  was  used  on  the  stem.  The  copper  block 
required  two  layers  since  the  temperatures  were  nearly  double  those  of  the  stem. 


The  convective  heat  transfer  coefficient  (h)  at  the  outer  surface  of  the  copper  block  due  to  natural 
convection  was  assumed  in  the  range  of  1-6  W/m  K  (Incropera,  2002).  The  thermal  conductivity 
of  the  insulation  is  0.04  W/mK  giving  a  critical  radius  for  the  insulation  of  0.007-0.04  meters 
(Appendix  A).  Figure  13  shows  the  insulation  and  convective  shield  used  to  reduce  the  heat  loss. 


Figure  13.  Insulated  Copper  Block  and  Convection  Shield  Surrounding  the  Test  Stand 
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3.1.7  Water  Chiller 


The  pump/chiller  combination  (Figure  14a,  Tek-Temp,  Model  TKD-200,  208  VAC  single  phase) 
supplies  distilled  water  to  the  test  stand.  The  pump  can  supply  a  source  pressure  up  to  100  psi. 
The  chilleT's  temperature  control  is  dynamic  creating  a  periodic  fluctuation  (which  could  be 
ignored  in  time-averaged  calculations).  The  chilleT's  capacity  was  4.4  kW  at  20  °C  (15,000 
BTU/hr). 


During  initial  testing,  it  was  found  that  the  chiller  temperature  and  pressure  supply  were  erratic. 
Coiled  nylon  tubing  (35  feet  x  1/2  in  internal  diameter)  was  added  (system  volume  increase) 
externally  to  dampen  the  fluctuations  in  water  temperature  and  pressure.  The  valve  used  for 
pressure  and  flow  adjustments  are  shown  in  Figure  14b. 


(a) 


(b) 

Figure  14.  Chiller/Pump  (a)  for  Cooling  and  Flow  with  Valve  (b)  for  Pressure  Adjustment 
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Other  notable  fluid  eomponents  were  a  filter  (Swagelok,  SS-6TF-40)  with  a  40  micron  filter 
located  up  stream  of  the  system.  To  adjust  the  flow,  two  valves  were  located  in  parallel  after  the 
flow  meter.  One  was  a  small  adjustment  valve  (Nuclear  Products  Company)  while  the  other 
(Nupro,  SS-4D4S)  offered  large  flow  adjustment.  The  two  valves  (Figure  15)  provided  22  g/s  of 
flow  at  50  psi  (inlet  pressure). 


Figure  15.  Two  Valves  for  Gross  and  Fine  Flow  Adjustment 
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3.2  Instrumentation 


Characterization  of  the  flow  through  the  heat  exchanger  and  the  heat  flux  required  measurement 
of  the  pressure  difference  across  the  heat  exchanger,  temperatures  along  the  copper  stem, 
temperatures  at  the  inlet  and  exit  of  the  heat  exchanger,  water  mass  flow  rate,  load-force  values, 
heater  voltage,  shunt  voltage,  and  heater  current  values.  In  addition  to  the  above  requirements, 
additional  thermocouples  were  attached  to  the  heater  body  and  located  in  the  chamber  for  an 
energy  balance.  A  HP  Agilent  34970A  DAQ  unit  receives  and  relays  the  signals  from  the 
sensors.  The  Agilent  unit  is  capable  of  monitoring  four  DAQ  boards  with  22  channels  per  board. 
For  this  experiment,  only  one  board  was  used.  LabView  software  (version  8.6.1)  was  used  to 
record  measured  values.  All  runs  were  set  to  have  a  one  second  delay  after  each  channel  scan; 
data  was  recorded  at  three  second  intervals.  Details  of  the  LabView  programs  are  given  in 
Appendix  F. 


3.2.1  Thermocouples 

There  were  thermocouples  located  on  the  copper  stem,  copper  block  body,  inlet  and  exit  of  the 
heat  exchanger,  and  within  the  test  stand  convective  heat  transfer  shield.  The  stem 
thermocouples  were  used  for  determining  the  heat  flux  transferred  to  the  heat  exchanger.  The 
encapsulated  T-type  thermocouples  (1/16  in  diameter)  located  at  the  inlet  and  exit  of  the  heat 
exchanger  were  used  to  measure  the  bulk  temperatures  to  calculate  the  enthalpy  rise  of  the  water. 
The  thermocouples  were  held  in  position  using  nylon  compression  fittings.  A  four-way  fitting 
locates  them  in  the  fluid  stream  (Figure  16). 


Figure  16.  Inlet  (or  Exit)  Flow  Thermocouple  with  Pressure  Port  in  the  Bottom  of  the  Four- 

Way  Fitting 
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The  six  T-type  thermocouples  which  passed  through  insulation  were  attached  to  the  copper  stem 
(Omega,  0.005  in,  butt  welded)  (Figure  17).  The  thermocouple  beads  were  secured  to  the  stem 
by  epoxy  (Omegabond,  1500  “F)  in  six  0.020  in  diameter  x  0.030  in  depth  holes  (Figure  18). 
Various  methods  were  tried  for  thermocouple  attachment.  Heating  the  copper  block  while  curing 
the  epoxy  and  vibration  using  a  sonic  bath  to  remove  air  pockets  were  both  tried.  The  best 
method  was  to  let  the  epoxy  air  dry  but  use  more  liquid  bonding  agent  than  suggested  by  the 
manufacturer.  Thermocouple  stress  relief  supports  (aluminum)  were  added  to  help  prevent 
breaking  of  the  thermocouples  while  moving,  calibrating,  and  testing  (Figure  19). 


Figure  17.  Thermocouples  (Circled)  Attached  to  the  Sides  of  the  Stem 
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X 


Figure  18.  Location  and  Position  References  for  the  Stem  Thermocouples 


Figure  19.  Stress  Relief  Clamps  for  Heater  Thermocouples 
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3.2.2  Load  Measurement 


The  first  load  apparatus  had  a  solid  mechanical  method  for  applying  a  load  to  the  heat  exchanger. 
This  arrangement  required  continuous  adjustment  during  warm-up.  The  apparatus  was  re¬ 
designed  to  incorporate  a  spring  for  warm-up  and  solid  contact  when  steady  conditions  were 
reached.  A  threaded  knob  at  the  top  of  the  apparatus  applied  a  load  through  the  upper  and  lower 
spring  seats  to  the  top  center  of  the  load  cell.  The  load  cell  set  on  the  plunger  transmitted  the  load 
to  the  Glas-Therm  block  (Figure  12).  The  load  cell  (Omega,  L8 100-200- 100)  had  a  maximum 
load  limit  of  100  Ibf  and  requires  10  volts  DC  excitation.  Figure  20  is  a  drawing  of  the  load  cell 
and  the  loading  apparatus.  The  load  cell  is  shown  in  position  between  the  lower  spring  seat  and 
the  plunger. 


Figure  20.  Load  Apparatus  and  Load  Cell  on  the  Upper  Alignment  Plate 


3.2.3  Pressure  Transducers 

The  pressure  transducers  (Omega,  PX303-050A5V)  had  a  50  psi  maximum  limit  and  required 
excitation  voltages  in  the  range  9-30  VDC  with  output  voltages  in  the  range  0.5-5. 5  VDC.  Figure 
21  shows  an  image  of  one  pressure  transducer  for  the  exit  flow.  Calibration  methods  are 
described  in  Appendix  C.  The  pressure  transducer  of  Figure  21  was  connected  to  the  system 
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through  a  4-way  Tee  at  the  exit  side  of  the  heat  exchanger  (flow  taps).  A  second  pressure 
transducer  was  located  behind  the  test  stand  of  Figure  21. 


Figure  21.  Pressure  Transducer  in  Front  of  the  Test  Stand 


3.2.4  Mass  Flow  Rate 

The  mass  flow  rate  of  the  water/cooling  fluid  was  measured  with  a  turbine  flow  meter.  The  flow 
meter  provided  a  three-wire  analog  output  to  a  signal  transmitter  (Omega,  FLSC61)  (Figure  22). 
The  flow  meter  uses  a  magnetic  wheel  as  a  signal  interrupter  to  produce  a  frequency  signal.  The 
signal  was  relayed  to  the  transmitter  where  it  was  converted  into  a  corresponding  output  voltage 
(millivolts).  Details  of  the  mass  flow  meter  calibration  method  are  provided  in  Appendix  C. 
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Figure  22.  Flow  Meter  with  Transmitter 


3.2.5  Electrical  Supply  Components 

The  power  supply  for  the  pressure  transducers,  flow  meter,  and  load  cell  was  provided  by  one 
unit,  a  Hyperion  Mode  DC  power  supply  (HY-WI-20-1.5)  (Figure  24).  The  pressure  transducers 
accepted  a  supply  voltage  in  the  9-30  VDC  range,  the  flow  meter  was  in  the  range  of  12-35 
VDC,  but  the  load  cell  was  specific  at  10  VDC.  The  power  supply  is  set  to  18  VDC.  A  voltage 
divider  is  used  to  drop  the  voltage  to  the  required  10  volts  for  the  load  cells  (Figure  23). 


Figure  23.  Schematic  of  the  Voltage  Divider  Installed  for  the  Load  Cell  Circuit 
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Figure  24.  Power  Supply  for  Pressure  Transducer,  Flow  Meter,  and  Load  Cell 


The  heater  power  supply  (Sorenson,  DCR  600-1. 5B)  is  800  VDC  and  2  A  capable.  Figure  25 
shows  the  power  supply,  and  Figure  26  shows  the  cylindrical  heater.  The  heater  (3/8  in  diameter 
X  2  in  length)  could  provide  200  W  of  power  at  120  volts.  The  heater  internal  resistance  was  76 
Ohms  and  requires  a  housing  ground.  Early  tests  showed  large  power  fluctuations  over  short 
time  steps.  To  help  with  this  issue,  two  capacitors  (2500  pF)  were  placed  in  line  to  smooth  the 
power  variation  (Figure  27). 


Figure  25.  Power  Supply  for  the  Cylindrical  Heater  within  the  Copper  Heater  Block 
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Figure  26.  Cylindrical  Heater  Located  within  the  Copper  Block  (on  the  Right) 


Power  Supply  . 

^  Capacitors  Heater 


Figure  27.  Schematic  of  the  Capacitors  Connected  to  the  Heater 
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4.0  EXPERIMENTAL  PROCEDURE 


The  procedure  to  run  an  experiment  included  adjusting  the  chiller  temperature  and  pressure,  flow 
rate,  contact  load,  and  heater  and  sensor  power.  Prior  to  performing  a  test,  the  electrical 
equipment  and  the  cooling/pump  were  run  20-30  minutes  prior  to  data  collection  to  allow  for  the 
temperature  of  the  involved  circuits  to  stabilize. 


4.1  Water  Chiller 

The  inlet  fluid  temperature,  20  °C,  was  set  at  the  chiller/pump  unit  using  the  available  PID 
(Proportional  Integral  Derivative)  controller.  The  fluid  temperature  measured  at  the  inlet  to  the 
heat  exchanger  varied  by  a  few  degrees,  depending  on  the  flow  rate  of  the  system.  The  inlet/exit 
heat  exchanger  thermocouples  (TC  13  &  14  -  channels  113  and  1 14)  and  the  number  1  stem 
thermocouple  (TC  1  -  channel  101)  were  monitored  for  steady  periodic  conditions.  The  steady 
periodic  conditions  were  selected  by  viewing  the  temperature  fluctuations  over  a  range  of  time 
(on  the  order  of  an  hour.  Figures  28  and  29).  The  steady  periodic  peak-to-peak  variation  ranges 
from  0.5°C  at  low  flow  (highest  AT  in  fluid)  to  below  0.1  °C  at  higher  flow  rates.  The  effect  that 
the  fluctuations  have  on  the  calculations  and  analysis  is  minimal  for  long  periods.  Figure  30 
shows  that  the  temperature  fluctuation  between  the  inlet  and  exit  of  the  heat  exchanger  is  small 
(0.04  “C  standard  deviation).  Also,  the  variation  of  the  inlet  fluid  temperature  would  have  a  small 
effect  on  the  stem  calculations  for  energy  flux  with  temperatures  of  100°C  at  TCI  and  VO^C  at 
TC3  (Figure  31). 
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Figure  28.  Sample  Temperature  Plot  for  Inlet  Fluid  Thermocouple 
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Figure  29.  Sample  Stem  Temperature  Trend  (Thermocouple  #1,  Same  Test  as  in  Figure  28) 
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Figure  30,  Temperature  Difference  between  Inlet  and  Exit  of  Heat  Exchanger 
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4.2  Pressure  and  Flow  Rate 


The  inlet  pressure  in  the  eoolant  lines  was  set  to  approximately  50  psi  (52  psi  max)  at  the 
ehiller/pump.  The  pressure  setting  was  adjusted  by  the  ball  valve  loeated  at  the  rear  of  the  unit. 
This  ball  valve  is  a  bypass  valve  whieh  diverts  unneeded  output  flow  baek  to  the  eooler/pump. 
With  the  system  set  at  50  psi,  a  maximum  flow  rate  of  22  g/s  was  possible. 


The  heat  exehanger  flow  rates  were  in  the  range  of  2-22g/s.  The  desired  flow  rates  were  obtained 
by  adjusting  the  low  flow  and  high  flow  needle  valves  loeated  downstream  of  the  flow  meter  and 
simultaneously  adjusting  the  line  pressure  (Figure  15).  A  large  needle  valve  is  used  for  gross 
adjustment  and  a  miero  valve  is  used  for  fine  adjustments.  Mass  flow  rates  below  3-3.5  g/s 
resulted  in  high  flow  fluetuations  and  were  diffieult  to  quantify. 


4.3  Contact  Load 

The  contact  load  (pressure)  at  the  interface  of  the  heat  exchanger  and  the  stem  was  relatively  low 
(20  lbf~  200  psi).  repeatable  conditions  were  important  for  testing.  This  load  is  suspected  to  be 
part  of  the  repeatability  problem.  It  was  later  found  that  the  load  was  insufficient  to  overcome  the 
friction  between  the  load  plunger  and  the  upper  plate  (Figure  12).  As  will  be  described  later,  the 
repeatability  problems  were  resolved. 


4.4  Power  Supplies 

The  electrical  power  applied  to  the  heater  within  the  copper  block  was  set  to  the  desired  rating  by 
adjusting  the  current  dials,  course  and  fine,  on  the  power  supply.  The  heater  resistance  was  76 
Ohms  at  room  temperature  and  increases  slightly  as  the  heater  warms.  The  voltage,  current,  and 
power  are  monitored  using  LabView  DAQ  software.  A  shunt,  installed  in  series  with  the  heater 
supply  power,  provides  a  voltage  drop  that  translates  to  a  current  value.  The  shunt  (0.02  Ohm) 
has  an  accuracy  of  0.31%  as  reported  by  the  manufacturer  (RC  Electronics).  The  voltage  drop  is 
measured  by  the  Agilent  DAQ  unit.  The  imposed  power  on  the  heater  drifts  with  time.  Thus, 
monitoring  and  adjustment  is  needed  periodically  during  runs.  The  only  requirement  by  the 
power  supply  for  the  load  cell,  pressure  transducers,  and  flow  meter  (Figure  24)  was  to  verify 
that  it  was  set  at  18  VDC  to  prevent  over-powering  the  load  cell. 
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5.0  RESULTS  AND  DISCUSSION 


A  test  stand  was  designed  to  provide  two  methods  for  energy  balanee  verification.  The  copper 
stem,  equipped  with  thermocouples  at  fixed  distances  along  two  sides,  provided  temperature 
values  used  to  determine  the  heat  flux  passed  to  the  heat  exchanger.  The  heat  flux  determination 
was  primarily  a  conduction  problem  and  is  considered  in  Section  5.1.  Analysis  of  the  convective 
heat  transfer  providing  the  enthalpy  rise  in  the  water  within  the  heat  exchanger  was  not  as 
straight  forward  as  the  conduction  heat  transfer  analysis  (Section  5.2).  The  convective  heat 
transfer  within  the  heat  exchanger  was  described  in  terms  of  Nusselt  numbers  and  comparisons 
are  made  with  established  correlations.  The  analysis  of  the  heat  transfer  includes  conduction 
along  the  nickel  heat  exchanger  material. 


5.1  Conduction  Heat  Transfer  through  Copper  Stem 

Two  calculations  associated  with  the  heat  exchanger  analysis  were  the  conductive  heat  flux 
through  the  stem  and  the  enthalpy  rise  in  the  working  fluid.  For  a  given  set  of  conditions, 
comparisons  of  the  enthalpy  rise  and  conductive  flux  provide  confirmation  of  the  heat  transfer. 
This  section  was  primarily  concerned  with  the  quantification  of  heat  conduction  through  the 
stem.  The  evaluation  of  the  heat  conduction  through  the  stem  was  performed  using  the  voltages 
(and  thus  the  temperatures)  provided  by  six  thermocouples  placed  on  opposite  surfaces  of  the 
stem  (Figure  31).  The  opposing  surfaces  correspond  to  the  inlet  and  exit  of  the  mini-channel  heat 
exchanger. 


Figure  31.  Six  Thermocouple  Positions  on  Copper  Stem  (Upper  Portion  of  Heater  Block) 


The  thermocouple  data  averaged  over  a  few  hundred  seconds,  after  reaching  a  steady-state 
condition  for  each  water  flow  rate  and  heater  power  input  setting,  were  analyzed  and  compared. 
Figure  32  is  a  plot  of  the  temperatures  obtained  from  the  thermocouples  located  at  the  relative 
position  number  on  the  stem  (Figure  3 1)  for  a  water  mass  flow  rate  of  7g/s  and  a  90  W  heater 
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input.  The  #1  and  #4  thermocouples  are  both  shown  as  position  1  in  the  plot  because  they  are  at 
the  same  level  on  stem  (same  distance  from  the  tip).  The  same  method  is  used  for  thermocouples 
2-5  and  3-6. 


Postion  # 

Figure  32.  Temperature  Distribution  along  Copper  Stem  for  Water  Flow  Rate  of  7g/s 


Initially,  it  was  thought  that  the  nonlinearity  of  the  profiles  of  Figure  32  were  due  to  uncertainty 
in  positioning  or  temperature.  However,  this  is  not  believed  to  be  the  case,  after  some 
investigation,  as  explained  in  the  next  several  paragraphs.  First,  the  influence  of  thermal 
expansion  on  the  final  location  the  thermocouples  were  determined  (i.e.,  differences  between  the 
two  temperature  profiles).  The  chemical  composition  of  the  stem  is  99.9%  oxygen  free  copper. 
Thus,  expansion  and  movement  of  a  thermocouple  may  be  possible  (Chemical  analysis  in 
Appendix  G).  Copper  has  a  thermal  expansion  coefficient  (a)  of  9.4  x  10'^in/°F.  The  stem  is 
approximately  8  mm  square  and  has  a  length  of  0.01 194  m  (L).  The  linear  expansion  due  to  an 
increase  in  temperature  is  given  by  (Holman,  2001): 

£  =  aATL.  (2) 

For  a  temperature  difference  of  200°C  (which  is  the  greatest  temperature  rise  expected),  the 
thermal  expansion  is  2.387  x  10'^  m.  This  value  has  a  negligible  influence  on  the  thermocouple 
location  and,  thus,  on  the  measured  axial  temperature  profde. 


Another  reason  for  a  nonlinear  profile  may  be  the  heat  loss  from  the  vertical  surfaces  of  the  stem. 
The  stem  contains  one  layer  of  insulation  with  a  k  value  (conduction  heat  transfer  coefficient)  of 
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0.04  W/m  K.  The  maximum  heat  transfer  eoefflcient  value  of  6  W/m^  °C  on  the  inside  of  the 
chamber  was  estimated  by  treating  the  stem  surfaces  as  flat  plates  in  the  presence  of  natural 
convection.  The  heat  flux  from  the  stem  sides  is  at  most  330  W/m  .  The  total  transfer  is  (stem 

3  2 

surface  area  =  4.064  x  10'  m  )  0.134  W.  Thus,  the  nonlinearity  in  the  stem  temperature  profile  is 
not  due  to  heat  loss  from  the  stem. 


Since  neither  the  heat  loss  from  the  stem  vertical  surfaces  nor  the  thermal  expansion  is 
responsible  for  the  measured  non-linear  temperature  profile,  another  analysis  was  performed  to 
consider  the  effect  of  thermocouple  placement  on  the  stem  temperature  profile.  There  could  be 
significant  error  due  to  machining  or  thermocouple  bead  placement.  Machining  error  often  is  a 
directional  error,  that  is,  if  a  hole  is  offset  in  one  direction  due  to  machining  error,  typically  all 
other  errors  will  be  offset  in  same  direction. 


Measurement  of  the  thermocouple  hole  positions  was  performed  after  machining.  An  average 
value  for  Ax  (0.0047  m)  is  used  for  the  following  analysis.  The  location  of  the  thermocouples 
within  the  0.020  in  diameter  holes  in  the  copper  stem  also  has  an  uncertainty.  The  bead  diameters 
of  the  thermocouples  (2  -  0.005  in  thermocouple  wires)  are  0.01-0.01 1  in  This  translates  to  a  + 
0.005  in  vertical  position  uncertainty  for  each  thermocouple  in  the  stem. 


Figure  33  shows  the  temperature  profile  along  the  stem  with  ±0.005  in.  thermocouple  placement 
indicated  by  horizontal  error  bars.  For  reference,  the  dashed  line  indicates  a  linear  temperature 
profile.  The  uncertainty  of  the  temperature  difference,  thermocouple  to  thermocouple,  is  less 
than  0.5  ‘’C.  The  uncertainty  is  the  standard  deviation  of  the  temperature  difference  for  the  fluid 
thermocouples  during  calibration  (no  statistical  analysis).  Figure  33  reveals  that  positioning 
uncertainty  for  the  thermocouples  would  not  correct  the  temperature  profile  along  the  copper. 
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Figure  33.  Stem  Temperature  Profile  (  +0.005  in  Thermocouple  Placement) 

The  temperature  profile  in  the  stem  is  better  understood  by  eonsidering  a  temperature  eontour 
plot  involving  the  heat  exehanger  obtained,  in  part,  by  CFD  ealeulations.  Figure  34  shows  a 
sketeh  of  the  expeeted  stem  temperature  profile  due  to  eontaet  with  the  heat  exchanger.  The 
temperature  distribution  within  the  heat  exchanger  (Figure  35)  was  calculated  by  Fluent,  a 
commercial  CFD  code.  The  stem  temperature  contours  were  visually  extrapolated  from  the  heat 
exchanger  temperature  contours. 


Fluid 


Figure  34.  Sketch  of  Expected  Heat  Conduction  Flux  Directions  within  the  Stem 
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This  contour  plot  suggests  that  the  heat  flux  calculation  on  the  stem  is  not  one-dimensional  near 
the  heat  exchanger.  The  high  heat  transfer  rate  at  the  upstream  contact  point  of  the  heat 
exchanger  and  the  stem  has  the  effect  of  causing  a  higher  conduction  heat  transfer  rate  at  that 
comer  of  the  stem.  Thus,  the  four  thermocouples  furthest  from  the  heat  exchanger  (1,  2,  4,  5  - 
Figure  31)  should  be  used  to  determine  the  heat  flux  through  the  stem.  The  conduction  rate  using 
the  thermocouples  is  72.48  W  (average  of  both  sides).  The  enthalpy  change  for  the  same 
condition  (7  g/s  flow)  is  72.66  W.  These  values  show  that  there  is  very  good  agreement  between 
the  energy  conducted  through  the  stem  and  the  change  in  enthalpy  of  the  cooling  fluid.  The 
energy  balance  for  the  stem  and  fluid  (2-22  g/s)  is  plotted  in  Figure  35.  The  solid  line  indicates  a 
1:1  balance.  The  heat  flux  passing  through  the  0.64  cm  (8x8  mm)  stem  at  the  7  g/s  conditions 
is  1 13  W/cm^.  This  meets  the  requirement  of  100  W/cm^.  In  pre-testing,  -160  W/cm^  heat 
transfer  was  observed  an  interface  temperature  (stem-heat  exchanger)  of  90  °C  at  7  g/s  flow  in 
the  heat  exchanger.  The  low  temperatures  at  the  interface  indicates  that  the  thermal  resistance 
between  the  stem  and  heat  exchanger  could  be  doubled,  from  the  insertion  of  SiC  device,  before 
the  temperature  limit  would  be  reached. 


Heat  Transfer  Through  Stem  (Watts) 


Figure  35.  Energy  Balance:  Stem  Conduction  and  Enthalpy  Rise  of  Water 
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5.2  Convective  Heat  Transfer  within  Mini- Channel  Heat  Exchanger 


In  Equation  6,  Cp  is  the  constant  pressure  specific  heat,  rh  is  the  mass  fiow  rate,  and  AT  is  the 
difference  between  the  bulk  inlet  and  outlet  temperatures. 

Figure  36  shows  the  calculated  enthalpy  change  in  the  heat  exchanger  fluid,  Q,  recorded  over 
three  consecutive  days  of  testing.  There  are  discontinuities  in  the  curve  at  7  and  1 1  g/s.  The 
discontinuities  are  believed  to  be  due  to  repeatability  issues  connected  with  the  contact 
resistance.  At  the  writing  of  this  report,  not  all  of  the  experiments  had  been  repeated  multiple 
times.  However  a  measure  of  confidence  in  the  values  comes  from  two  months  of  setup  and 
validation  tests.  Before  final  testing  began  (for  the  results  presented  here),  validation  tests  were 
run  to  help  with  test  stand  adjustments.  These  tests  were  all  performed  at  7  g/s.  The  heat  transfer 
values  and  temperature  conditions  found  in  the  first  group  of  data  (before  the  first  discontinuity) 
are  very  close  to  the  conditions  recorded  in  the  pre-tests.  Moreover,  the  maximum  difference  in 
the  heat  transfer  rates  at  the  discontinuities  of  1 1  g/s  is  below  2%.  At  flow  rates  of  20  g/s  and 
greater,  the  heat  transfer  rate  to  the  heat  exchanger  would  normally  be  expected  to  increase  when 
the  Rcd  are  in  the  range  of  laminar  to  turbulent  flow  transition.  However,  there  is  a  decrease  in 
heat  transfer.  This  decrease  in  the  heat  transfer  rate  near  20  g/s  is  believed  to  be  a  result  of  the 
plunger  (last  component  in  load  apparatus  transferring  force  to  the  heat  exchanger)  sticking.  The 
repeatability  issues  and  potential  plunger  sticking  problems  would  need  to  be  further  explored  if 
more  work  was  undertaken. 
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Figure  36.  Heat  Transfer  Rate  to  the  Eleven  Channels  of  the  Heat  Exchanger 


The  surface  temperature  for  each  channel  within  the  heat  exchanger  is  unknown,  so  an  estimate 
is  used.  Using  the  thermocouple  values  in  the  copper  stem  and  the  thermal  properties  for  the 
copper  and  nickel  (heat  exchanger),  a  temperature  was  estimated  at  the  bottom  of  the  channels. 
This  value  was  averaged  with  the  bulk  temperature  to  estimate  the  surface  temperature  for  the 
channels  and  assumed  the  same  temperature  around  each  channel.  Figure  37  shows  a  two- 
dimensional  color  contour  plot  for  a  plane  along  a  channel  axis  that  was  obtained  from  a  three- 
dimensional  CFD  simulation  (Fluent  code).  This  simulation  assists  in  the  explanation  of  the 
convective  heat  transfer  coefficient  calculation  and  is  useful  for  comparison  of  the  results  with 
existing  heat  transfer  correlations.  The  plane  for  the  color  contour  plot  shown  is  located  at  the 
mid-point  of  the  heated  section.  Figure  37  shows  that  the  surface  temperature  varies  relatively 
little  from  3 1 .4  “C  to  32.2  “C  with  an  average  of  3 1 .8  ±  0.4  °C. 
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Figure  37.  Calculated  Temperatures  (K)  inside  and  outside  the  Center  Channel  (Rcd  of 

1260) 


The  value  for  L  (heat  transfer  length  of  the  channels  in  Eqn.  4)  is  important  for  accurate  results. 
As  seen  in  Figure  34  the  heat  exchanger  extends  beyond  the  stem  and  acts  to  spread  the  heat  by 
conduction.  The  heat  transfer  length  selected  to  represent  the  region  of  convective  heat  transfer 
could  be  the  channel  length  (gives  a  low  Nud  value),  stem-to-heat  exchanger  contact  length,  or 
more  appropriately  the  actual  length  along  the  channel  in  which  the  heat  is  conducted  in  the  solid 
wall  (i.e.,  conjugate  heat  transfer  effects  obtained  by  CFD  simulation). 


Figure  38  is  a  plot  of  Nud  with  Rcd  for  two  of  the  above  situations  -  shortest  and  longest  heat 
transfer  lengths  (8  mm  and  21.59  mm).  Nud  is  defined  as  (Incropera,  2002) 


Nuo  =  -Y  (7) 

In  Fqn.  7,  D  is  the  diameter  of  a  channel,  and  k  is  the  thermal  conductivity  of  the  water.  The 
convective  heat  transfer  coefficient,  h,  is  determined  from  Fqn.  3.  The  Reynolds  number  was 
based  on  average  velocity  at  the  inlet  temperature  (which  varies  relatively  little  from  the  bulk 
temperature)  and  is  defined  as 


Rbd 


pDV 


(8) 


where  p  is  the  fluid  density  ,  the  diameter  of  the  channels,  V  is  the  average  velocity,  and  p  is  the 
dynamic  viscosity  of  water. 
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Figure  38.  Nud  Resulting  from  Use  of  the  Shortest  and  Longest  Heat  Transfer  Lengths 

Figure  38  shows  that  a  Nud  profile  for  Rod  in  the  range  0-1900  is  followed  by  a  step  ehange  in 
the  profile  for  larger  Ren  (a  step  at  1100  Re  exists  but  is  mueh  smaller).  The  raw  data  for  the 
three  segments  of  a  Nud  curve  were  obtained  on  consecutive  days.  The  sudden  change  in  a  Nud 
curve  is  due  to  repeatability  issues  which  vary  with  daily  setup.  The  decrease  in  all  of  the  Nud 
curves  near  Rcd  of  3800  is  believed  to  be  due  to  load  plunger  issues.  These  repeatability  issues 
were  described  earlier  in  connection  with  Figure  36  correspond  to  a  decrease  of  approximately 
13%  and  are  not  considered  fatal  to  interpretation  of  the  experimental  results. 


CFD  simulations  were  used  to  determine  more  appropriate  values  for  the  heat  transfer  length.  A 
two-dimensional  model  was  used  to  represent  the  maximum  thermal  spreading,  and  a  three- 
dimensional  model  was  used  for  a  mid-range  flow  rate  (7  g/s).  The  minimum  heat  transfer  length 
is  the  length  that  the  stem  makes  contact  with  the  heat  exchanger  (8  mm).  Figure  39  shows  the 
contact  interface  of  the  stem  and  the  heat  exchanger.  The  stem  is  slightly  recessed  into  the  heat 
exchanger.  The  heat  exchanger  and  stem  together  form  a  line  of  contact  parallel  to  the  fluid  flow. 
This  line  of  contact  is  8mm  in  length. 
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Figure  39.  The  Heat  Exchanger  and  Stem  Form  a  Line  of  Contact  Parallel  to  Fluid  Flow 


Figure  40  shows  a  contour  plot  of  the  channel  region  calculated  by  the  commercially  available 
CFD  code,  CFDACE.  For  a  flow  rate  of  2  g/s  and  a  heat  flux  of  70  W/cm  ,  the  two-dimensional 
numerical  simulation  shows  how  the  channel  wall  thickness  allows  the  imposed  heat  flux  to 
diffuse  such  that  the  heated  length  is  not  the  same  as  the  contact  length  between  the  stem  and 
heat  exchanger.  The  estimated  length  from  the  CFD  simulation  is  ~  14. 5  mm  and  will  vary  with 
the  heat  transfer  and  mass  flow  rate  in  a  channel.  To  determine  the  heated  characteristic  length 
from  the  temperature  contour  plot,  an  examination  of  the  thermal  boundary  layer  is  required. 
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Figure  40.  Temperature  Contour  Plot  for  a  Two-Dimensional  CFD  Solution  (Rcd  of  1260) 
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The  thermal  boundary  layer  is  considered  here  to  be  the  region  where  there  are  significant 
differences  (more  than  3  °C)  between  the  fluid  and  channel  surface  temperatures.  The  length  of 
the  thermal  boundary  layer  was  taken  as  L.  A  three-dimensional  simulation  was  also  performed. 
For  the  three-dimensional  length  estimate,  axial  and  radial  temperatures  were  taken  into  account. 
For  a  Rcd  of  1260  (mass  flow  rate  of  7  g/s),  the  measured  bulk  temperature  increase  of  the  fluid 
from  the  inlet  to  exit  of  the  heat  exchanger  was  2.46°C.  The  three-dimensional  model  gave  a 
temperature  increase  of  2.5°C  for  the  same  conditions  which  is  in  very  good  agreement  with  the 
measurements.  Unfortunately,  the  three-dimensional  model  required  a  large  numbers  of  cells  and 
extended  run  times  were  required.  For  this  reason,  the  three-dimensional  simulations  were 
performed  only  at  7  g/s  (1260  Rcd).  The  analysis  for  the  heated  length  for  this  model  provides  a 
length  of  approximately  1 1mm.  Figure  41  shows  the  temperature  contour  plot  for  the  conditions 
tested. 
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Figure  41.  Two-Dimensional  Temperature  Contours  (from  Three-Dimensional  Model) 

along  the  Channel  Centerline  (Reo  of  1260) 


Not  having  multiple  contour  plots  meant  that  the  temperature  profiles  for  all  the  experimental 
flow  rates  would  not  be  available.  To  fill  the  gap  for  all  other  flow  rates  tested,  an  exponential 
profile  is  assumed  to  fit  the  lengths  as  a  function  of  flow  rate  for  simplicity.  The  exponential 
curve  was  fit  to  three  values:  (the  shortest  expected  length)  8mm  (22g/s),  1 1mm  (7g/s  from  a 
three-dimensional  CFD  model),  andl4.5mm  (2g/s,  the  lowest  flow  rate  from  a  two-dimensional 
CFD  model).  The  values  for  the  heated  lengths  were  then  estimated  for  each  of  the  flow  rates. 


Figure  42  is  a  plot  of  Nud  with  Rcd  which  is  determined  from  the  average  velocity  for  water  in 
the  channels  and  the  mean  bulk  temperatures.  The  Nud  number  is  determined  with  an  area  based 
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on  a  variable  heat  transfer  length  (eurve  fit).  The  Nud  with  variable  heat  transfer  length  is 
compared  to  Nud  for  constant  lengths  corresponding  to  the  channel  length  and  the  contacting 
interface  length  (between  the  heat  exchanger  and  copper  stem).  The  plot  indicates  that  heat 
diffusion  in  the  walls  of  channels  and  tubes  may  strongly  affect  Nud  and  may  need  to  be  taken 
into  account.  For  the  situation  of  a  constant  heat  flux  with  fully  developed  laminar  flow,  Nud  is 
4.36.  We  do  not  have  fully  developed  flow  and,  thus,  NuD  should  be  above  4.36  as  in  Figure  42. 


Reo 

Figure  42.  Comparison  of  Three  Characteristic  Heat  Transfer  Lengths 


Comparisons  with  established  heat  transfer  and  flow  correlations  assist  in  validation  of  the 
assumptions  made  with  respect  to  the  heat  transfer  length  described  in  preceding  paragraphs. 
Heat  transfer  correlations  are  available  for  various  entry  conditions:  combined  thermal  and 
hydrodynamic,  only  hydrodynamic,  or  only  thermal.  Each  of  the  conditions  is  available  with 
either  uniform  heat  flux  or  constant  temperature  at  the  channel  wall.  Selected  correlations  are 
now  considered.  The  hydrodynamic  boundary  length  is  calculated  as  (Wilcox,  2000) 

^  =  Q.QSReD  (9) 

The  X  value  (fully  developed  length)  is  11.3  mm  at  the  lowest  laminar  mass  flow  rate  to  1 14  mm 
at  the  highest  (358  <  ReD  <  3600).  For  turbulent  flow  (3600  <  Rcd  <3900),  the  length  of  the 
entry  region  is  described  by  (Wilcox,  2000) 
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(10) 


The  length  (x)  associated  with  a  diameter  of  0.000635  m  is  0.00635  m  or  roughly  1/3  of  the 
channel  length.  It  is  important  to  know  where  the  temperature  profile  will  become  fully 
developed.  An  estimate  of  the  thermal  entry  length  is  given  by  (Wilcox,  2000) 

^-O-OSKeoPr  (11) 

In  Eqn.  II,  the  Prandtl  number  (6.996)  is  evaluated  at  the  mean  bulk  temperature.  The  estimated 
entry  lengths  are  0.0795  to  0.866  m. 


5.2.2  Comparison  with  Experimental  Data 

The  first  comparison  between  the  experimental  data  and  a  heat  transfer  correlation  is  for  flow 
with  developing  thermal  and  hydrodynamic  profiles  (Eigure  43).  A  correlation  for  the  mean 
Nusselt  number  (Shah  (1978))  used  for  small  channels  with  a  combined  entry  length  and  uniform 
heat  flux  is  (Young,  2008) 


Nuo=l  ^  .  (12) 

f  or  Re  Pr(—j  >  33.3 

To  illustrate  the  consequences  of  the  selection  of  E,  E  in  Eigure  43  is  chosen  as  the  entire 
channel  length  (21.59  mm)  for  simplicity  rather  than  the  heat  transfer  length  derived  from  CED 
simulations.  Thus,  the  entire  channel  length  is  selected  for  both  the  evaluation  of  Eqn.  12  and  the 
experimentally  determined  Nud.  Here,  Rcf  is  the  Reynolds  number  based  on  an  average  fluid 
velocity,  and  the  Prandtl  number  is  obtained  using  the  mean  bulk  temperature.  The  bulk  viscosity 
(from  the  mean  bulk  temperature)  and  the  estimated  viscosity  at  the  wall  are  |Uf  and  |Uw, 
respectively.  With  this  choice  of  heat  transfer  length,  Eqn.  12  appears  to  overestimate  Nud  values 
relative  to  those  determined  experimentally.  We  know  from  the  CED  simulations  that  the  entire 
channel  length  is  not  the  optimum  choice  for  E.  Thus  without  the  CED  simulation,  it  is  not  clear 
a  priori  what  value  of  E  should  be  used  in  the  experimental  data  reduction  or  for  use  with  heat 
transfer  correlations. 
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Figure  43.  Nud  (Shah  (1978))  Correlation  for  Combined  Thermal  and  Hydrodynamic 
Entrance  Length  with  Nud  Determined  from  Experimental  Measurements  (L  =  21.59  mm) 

Comparison  of  Nud  that  are  experimentally  obtained  here  with  those  from  heat  transfer 
eorrelations  available  in  the  literature  ean  provide  additional  eonfirmation  of  the  measurement 
methodology.  That  is,  provided  L  from  the  CFD  caleulations  (or  a  more  judicious  choice  than 
simply  guessing)  is  chosen  for  use  in  reduction  of  the  measured  data  when  its  value  is  otherwise 
unclear.  Most  available  heat  transfer  correlations  for  the  flow  and  thermal  entry  region  within  a 
circular  tube  were  derived  for  either  a  uniform  heat  flux  or  isothermal  wall  boundary  along  the 
entire  tube.  This  is  not  the  case  for  the  present  experiments  as  the  tube  is  heated  only  over  a 
section  of  the  tube  and  only  over  a  portion  in  the  circumferential  direction.  In  addition  with 
regard  to  the  entry  region,  heat  transfer  correlations  are  available  for  developing  temperature 
profiles  assuming  a  fully  developed  velocity  profile  and  for  velocity  and  temperature  profiles 
which  simultaneously  develop  from  the  inlet.  In  this  work,  the  heated  tube  section  is  downstream 
from  the  inlet,  and  Pr  is  ~7.  Thus,  the  thermal  and  velocity  profiles  do  not  develop  at  the  same 
rate,  and  the  starting  location  of  the  thermal  boundary  layer  does  not  coincide  with  that  of  the 
velocity  boundary  layer.  Although  the  requirements  of  the  available  heat  transfer  correlations  do 
not  quite  match  the  actual  experiments,  it  is  expected  that  they  will  still  bracket  the  experimental 
Nud  values. 


For  conditions  of  a  uniform  heat  flux  and  for  thermal  and  velocity  profiles  which  develop  from 
the  inlet,  a  local  Nud-x  correlation  (local  axial  location  denoted  by  the  subscript  -x)  is  given  by 
Churchill  and  Ozoe  (Bejan,  1995)  (Churchill,  1973): 
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Here,  the  Graetz  number  (Gz)  is  defined  as 


TT 


(14) 


and  the  Prandtl  number  for  the  mean  bulk  fluid  temperature  (Pr)  is  6.996.  x*  is  defined  as 


jx/D) 

Rcj)  Pr 


(15) 


Equation  13  was  numerieally  integrated  from  the  channel  entrance  to  the  end  of  heated  copper 
contact  region  (from  x  =  0  mm  to  x  =  14.795  mm)  to  obtain  average  values  of  the  Nusselt 

number  (i.e.,  Nu^  ^ve  “  7  Jq  Nu0_xdx  )  for  a  range  of  Rcd.  Although  not  entirely  physically 

correct,  it  is  reasonable  to  choose  L  =  14.795  mm  for  use  with  Eqn.  13  for  simplicity.  Eor 
purposes  of  comparison,  the  Nud  values  obtained  from  the  integration  of  Eqn.  13  are  plotted 
together  with  the  experimental  Nud  values  (which  also  use  E  determined  from  CED  simulations 
and  referred  to  as  a  variable  heat  transfer  length  in  the  legend)  in  Eigure  44.  The  two  curves 
compare  reasonably  well  considering  the  typical  uncertainty  in  a  heat  transfer  correlation  is  + 
25%. 
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Figure  44.  Nud  Derived  from  the  Churchill  and  Ozoe  Correlation  (1973)  Plotted  with 

Measured  Values 


For  a  thermal  entry  region  alone,  assuming  a  fully  developed  velocity  profile  (Flagen-Poiseuille 
flow),  a  correlation  for  Nud-x  for  uniform  heat  flux  is  (Bejan,  1995;  Shah  and  London,  1978), 


NUd-x  =  < 


3.302x*-i/^  -  1.00  forx*  <  0.00005 
1.302x*-i/^  -  0.5,  for  0.00005  <  x*  <  0.0015 
4.364  +  8.68(10^x*  exp(-41x*)  ,  for  x*  >  0.001  ■ 

xfD 


(16) 


x*  =  ( 


RonPr 


;) 


Eqn.  16  is  essentially  a  fit  of  the  classic  Graetz  solution  (Kakac  et  ah,  1987).  Eqn.  16  was 
numerically  integrated  along  the  heated  copper  contact  length  (from  x  =  6.795  mm  to  x  =  14.795 
mm  with  the  inlet  as  the  origin)  to  yield  Nud  for  a  range  of  Rcd-  The  integration  range  used  for 
this  correlation  is  different  from  that  used  for  Eqn.  12  because  of  the  assumptions  of  fully 
developed  flow  (which  does  not  actually  exist)  and  a  developing  thermal  profde  (along  this 
length).  This  choice  of  the  integration  limits  emphasizes  the  effect  of  the  developing  thermal 
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profile.  (In  Eqn.  16,  Pr  and  Rcd  are  evaluated  at  the  mean  bulk  temperature.)  The  Nud  values 
obtained  from  the  integration  of  Eqn.  16  and  the  experiments  are  plotted  together  in  Eigure  45. 
Eigure  45  shows  that  the  integration  of  this  correlation  over  the  selected  integration  limits 
underestimates  Nud  obtained  from  the  experiments  (which  use  E  from  the  CFD  simulations) 
over  this  Rcd  range. 


Figure  45.  Comparison  of  Shah  and  London  (1978)  Uniform  Heat  Flux  Correlation  with 

Experimental  Data 


A  heat  transfer  correlation  (Kakac  et  ah,  1987)  for  the  mean  Nusselt  number,  Nud,  that  was 
originally  derived  for  flat  ducts  has  also  been  used  for  circular  ducts  (Bejan,  1995)  is  given  by 


Nud  =  7.55  + 


0.024  X* 

(1  +  0.0358  PrO  i^x* -0-64)  ■ 


(17) 


This  correlation  is  for  simultaneously  developing  temperature  and  velocity  profiles  within  a  tube 
with  isothermal  wall  conditions,  and  x*  and  Pr  are  defined  in  the  same  way  as  in  Eqn.  16.  Eqn. 
17  was  evaluated  at  x  =14.795  mm  (the  end  of  the  heated  copper  contact  region)  to  provide  Nud 
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for  comparison  with  those  obtained  from  experiments  for  different  Reo  in  Figure  46.  The 
correlation  for  isothermal  walls  (Eqn.  17)  agrees  well  with  Nud  derived  from  the  present 
measurements. 


Figure  46.  Comparison  of  Experimental  Data  with  a  Correlation  for  Entry  Flow  with 
Constant  Surface  Temperature  (Kakac  et  al.,  1987) 

To  summarize,  Figure  47  shows  Eqns.  12  to  17  together  with  the  experimentally  acquired  Nud. 
Eqns.  12  and  17  are  the  best  fits  to  the  experimentally  determined  Nud  (variable  heat  transfer 
length).  However,  all  the  eorrelations  compare  reasonably  well  with  the  experimental  values. 
This  agreement  among  Nud  obtained  by  different  methods  provides  a  measure  of  confidence  in 
the  results.  (It  is  mentioned  in  passing  that  for  Reo  values  above  700  to  those  in  the  turbulent 
range,  the  fully  developed  hydrodynamie  region  is  beyond  the  channel  length.)  The  measured 
heat  transfer  values  involved  physical/mechanical  errors  which  resulted  in  values  7-11%  lower 
than  expected.  The  contact  resistance  for  the  thermal  grease  was  nearly  5  times  higher  than  the 
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manufacturer  specified,  indieating  poor  eontact.  This  is  known  because  losses  in  the  stem  are 
negligible.  The  uncertainty  in  Nud  is  addressed  in  Appendix  B.  However,  statistical  aspects  of 
the  uncertainty  are  not  fully  accounted  for  due  to  limited  testing  repetitions  (Note:  the  first  part 
of  the  experimental  data  is  used  for  all  eomparisons.) 
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6.0  FLOW  CHARACTERIZATION 


6.1  Head  Loss 

The  pressure  drop  across  the  heat  exchanger  was  measured  for  all  flow  passages  between  the 
inlet  and  exit  pressure  sensors  of  Figure  1.  It  is  important  to  understand  the  influence  of  each 
flow  element  on  the  available  pressure.  The  heat  exchanger  channels  themselves  are  of  interest 
for  fluid  flow  characterization  but  are  located  between  the  inlet  and  exit  pressure  transducers 
with  other  attachments  or  accessories.  The  measurement  of  pressure  drop  (AP)  across  a  channel 
section  can  be  used  to  evaluate  the  flow  losses.  This  means  that  the  AP  measurements  obtained 
here  will  necessarily  be  greater  than  those  of  the  heat  exchanger  channels  alone.  To  compare  the 
experimentally  determined  friction  factors  of  the  heat  exchanger  tubes  with  existing  friction 
factor  correlations,  for  example,  the  flow  effects  due  to  the  heat  exchanger  accessories  need  to  be 
removed.  Figure  48  shows  a  sectioned  view  of  the  heat  exchanger  (with  connecting  flow 
passages)  used  in  the  experiments.  Pressure  taps  are  located  at  either  end  of  the  (yellow)  tubes, 
and  the  letters  correspond  to  the  positions  listed  in  Table  1. 


Figure  48.  Sectional  View  of  Heat  Exchanger  Showing  the  Flow  Path 


Subtraction  of  the  pressure  drops  resulting  from  components  other  than  the  channels  is  now 
considered.  To  account  for  all  pressure  losses  in  the  system  Eqn.  18  is  used  (Maynes,  2002). 


^Ptotal  ~ 


KiV^ 


(18) 


Here,  p  is  the  density,  V  is  the  velocity, /is  the  friction  factor  for  a  channel  (a  function  of  tube 
surface  roughness  and  diameter),  L  is  the  component  length,  D  is  the  characteristic  diameter,  and 
Kl  is  the  pressure  loss  coefficient  which  accounts  for  losses  due  to  bends  and  expansions,  for 
example. 
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The  pressure  loss  coefficient  is  defined  as  Kl=  APl/(pV  /2)  where  AP  is  unknown  for  each 
section.  Kl  is  an  empirical  constant,  and  tabulated  values  for  different  flow  configurations  are 
available  (Idelchik,  1986).  Figure  49  shows  where  component  pressure  drops  occur  (Table  1). 


The  following  paragraphs  describe  the  path  through  the  channels  with  short  descriptions  on  how 
each  flow  situation  (change  in  flow  direction,  expansion,  contraction,  etc.).  At  the  ends  of  the 
two  tubes  (yellow)  of  Figure  48,  there  are  tees  with  ports  (not  shown)  that  are  connected  to 
pressure  transducers  by  tubing.  Starting  at  the  left  side  of  Figure  48,  the  flow  moves  through  the 
A”  tubing  to  the  mixing  chamber.  The  pressure  drop  in  this  short  tube  is  neglected.  At  this  point, 
the  flow  enters  the  inlet  of  the  mixing  chamber  (two  pieces  welded  -  blue/purple).  The  area  of 
the  inlet  tube  is  2.44  mm^  while  the  area  of  the  mixing  chamber  using  the  hydraulic  diameter 
(Dh)  is  4.4  mm  .  The  hydraulic  diameter  is  (Idelchik,  1986) 


4  *  ah 
2(a  +  b) 


(19) 


The  hydraulic  diameter,  the  equivalent  diameter  for  a  rectangular  channel,  is  a  function  of  the 
width,  a,  and  the  diameter,  b.  This  is  the  first  event  -  an  expansion  due  to  a  change  in  flow  area. 
After  the  flow  enters  a  chamber,  it  must  change  direction  (B  of  Figure  48).  This  event  (change  in 
direction)  can  be  represented  as  an  elbow.  The  elbow  flow  occurs  again  at  the  entrance  to  the 
channels.  Instead  of  two  separate  events,  the  flow  through  the  mixing  chamber  is  treated  as  a 
double  elbow,  one  event.  For  simplicity,  the  bend  radius  for  the  flow  is  estimated  at  14  the 
chamber  width.  This  changes  the  hydraulic  diameter  for  the  flow.  The  new  hydraulic  diameter  is 
2.54  mm^. 


The  entrance  to  the  heat  exchanger  channels  is  the  next  event.  This  location  has  two  flow  effects 
to  consider.  First,  there  is  a  contraction  from  the  mixing  chamber  to  the  heat  exchanger  channels. 
In  addition,  the  channels  have  a  sharp  entrance.  The  same  analysis  is  used  for  the  exit  side  of  the 
system  where  there  is  an  expansion  from  the  channels  to  the  mixing  chamber,  another  double 
elbow,  and  two  events  again  at  the  exit  of  the  mixing  chamber  to  the  A  tubing. 


Table  1  summarizes  the  flow  events  and  their  contribution  to  the  total  AP.  The  equations  in  the 
table  are  used  to  estimate  Kl  values.  With  the  Kl  coefficients,  the  head  loss  is  calculated  by 
dividing  Equation  18  by  p.  The  head  loss  for  each  flow  event  is  plotted  with  Rcd  within  the 
channels  in  Figure  49  using  the  average  velocity  which  was  calculated  using  transport  properties 
determined  at  20  “C.  Figure  49  shows  that  the  head  losses  are  greatest  for  the  entrance  and  exit 
regions  of  the  channels. 
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Figure  49.  Head  Loss  Values  (Idelchik,  1986)  Plotted  with  Reo  for  Each  Flow  Situation 


54 


Table  1.  Equations  (Idelchik,  1986)  Used  to  Calculate  Loss  Coefficients  for  a  Particular  Flow  Event  Neglecting  ^fr  (Minor  Loss 

Coefficient) 


Position  "Event"  Book  Page  K  (g) Components  of  Equations_ Details 


A 

"sudden  expansion  of  a  flow 
having  unifonn  velociw 
distribution" 

160 

(l-Fo'T:)'  +  ^fr'n' 

Fo  -  Inlet 
Area 

F2- 

Expansion 

Area 

^fr — Coeff. 
Of  friction 
resistanc 

e 

n=ratio  of 

areas 

F0=16.417mm^ 

F2=77.55mm^ 

n=3.76 

Vfr~0.017 

B 

"Z-sh^ed  elbows  with 
rounded  comers" 

306 

^  Savs  ^  loc  Sfr 

1 

C 

"entrance  into  a  straight 
tube  of  constant  cross 

section  with  tube  mounted 

flush  to  wail" 

122 

q(5/DK)  &  q(b/Dh) 

d  = 

thickness 

of 

extended 

tube 

b  =  length 
tube 
extends 
from  wall 

Dh  = 

hydraulic 

diameter 

d=0 

b=0 

Dh=635/im 

D 

"Z-shaped  elbows  with 
rounded  comers" 

306 

^  Savs  S  loc  Sfir 

E 

"sudden  expansion  of  a  flow' 
having  uniform  velocitv' 
distribution" 

160 

(I-Fq/F:)’  +  ^fr'n" 

Fo  -  Inlet 
Area 

Fz- 

Expansion 

Area 

F0=0.317mm^ 

F2=1.616mm^ 

n=3.76 

Vfr~0.017 

Dh=1  434mm 

F 

"sudden  contraction  in 

transition  and  laminar 
regions" 

168 

q(Fo/F  j.  Re)  graph 

FI  -  Inlet 
Area 

Fo- 

Contractio 
n  Area 

~=  -0.3151Ln(x)  + 2.6299 
(300-2000  Re) 

~=  0.3359Ln(x)-2  1565 
(20004000  Re) 

K  (i^)=AP  (pw'o*  2)  wo=veloctiy  used  for  head  loss  calculation  Cfr=Coeff.  Of  friction  -?''02  l^2h) 
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6.2  Characterization  of  the  Channels  -  Pressure  Loss 


Mini-  and  micro-  channels  provide  enhaneed  heat  removal  but  at  the  expense  of  increased  AP 
(pressure  differenee)  so  proper  optimization  must  be  considered.  Obtaining  accurate  pressure 
differenees  by  experiment  or  numerieal  simulation  is  therefore  eritieal  for  thermal  designs  if 
pumping  energy  savings  are  desired.  As  described  previously,  the  inlet  and  exit  pressure 
measurements  were  taken  across  the  entire  heat  exchanger  assembly,  whieh  includes  tube 
inlets/exits,  mixing  ehambers,  and  heat  exchanger  channels.  The  pressure  losses  due  to  the 
connecting  components  were  estimated  and  subtracted  from  the  AP  measured  from  the  test  stand. 
A  corrected  AP  for  the  heat  exehanger  channels  was  then  obtained  and  used  in  the  following 
analysis.  Although  there  are  uneertainties  in  this  indirect  method  (Appendix  2.6),  it  is  a  well- 
known  engineering  technique  and  reasonable  results  can  be  obtained.  Correlations  for  estimating 
the  departure  from  Hagen-Poiseuille  flow  (White,  1974)  and  the  turbulent  transition  flow  as  a 
function  of  Reynolds  number  were  used  to  help  validate  the  ealculated  friction  factor.  In  the 
current  work,  the  channel  flow  behavior  is  described  by  the  friction  factor,  /.  The  pressure  drop 
across  the  ehannels  can  be  used  to  determine /(Wilcox,  2000); 

2AP 


In  Eqn.  20,  V  is  the  loeal  average  veloeity.  All  Reynolds  numbers  (Reo)  are  based  on  the  average 
velocity  and  mean  bulk  temperature  of  the  water.  Figure  50  shows  the  measured  pressure  drop 
across  the  heat  exchanger,  the  measured  pressure  drop  minus  the  pressure  drop  from  head  loss 
calculations,  and  the  pressure  drop  given  by  the  Poiseuille  equation. 
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Figure  50.  Measured  AP  across  the  Heat  Exchanger,  Measured  AP  Minus  the  AP  from 
Head  Loss  Calculations,  and  the  AP  Given  hy  the  Poiseuille  Equation 


Evaluating  the  channel  friction  factor  (/)  for  this  case  has  a  degree  of  uncertainty  because  a 
measurement  of  the  surface  roughness  for  the  channels  was  unavailable.  The  channels  were 
manufactured  at  Oak  Ridge  National  Laboratory  using  an  electroplating  method  in  which  nickel 
was  accumulated  on  fiber  wires.  The  fibers  are  chemically  dissolved  after  the  process  is  finished. 
Figures  51a  and  b  are  scanning  electron  microscope  (SEM)  images  of  the  heat  exchanger 
channels.  The  images  show  that  the  channels  (635  |um  diameter)  have  smooth  circular  edges.  The 
range  of  surface  roughness  (s)  for  these  channels  is  approximated  as  15nm  to  1.5  |um  (a  smooth 
drawn  tube)  (Papautsky,  1999).  The  s/D  (ratio  of  surface  roughness  to  tube  diameter)  is  a 
maximum  of  0.0024  mm/mm.  Elsing  a  Moody  diagram,  the  expected  friction  factor  for  this  range 
is  ~0.038-0.042. 
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Figure  51.  SEM  Images  of  (a)  a  Single  and  (b)  Multiple  Channels 

The  calculated/ for  the  heat  exchanger  with  and  without  the  head  loss  due  to  the  flow 
accessories  is  shown  in  Figure  52.  The  adjustments  for  head  loss  place  the  calculated /near  the 
Hagen-Poiseuille  flow  line  (White,  1974),  as  is  should  be  for  a  smooth  channel.  One  concern  is 
the  first  section  of  the  experiment  data  line  (friction  factor  without  head  loss).  At  low  Rcd,  the 
curve  is  below  the  theoretical  curve.  This  is  believed  to  be  due  to  a  larger  uncertainty  in  the  mass 
flow  rate  at  low  Rcd  values.  However  since  the  channels  are  very  short,  the  inlet  and  exit  effects 
for  low  flow  rates  may  not  be  totally  accounted  for  with  the  present  head  loss  calculations. 
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Figure  52.  Comparison  of  Experimental  Friction  Factor  with  and  without  Head  Loss  from 
Accessory  Fittings  to  the  Darcy  Friction  Factor  for  Hagen-Poiseuille  Flow  (White,  1974) 
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A  comparison  was  made  between  the  current  experimental  data  and  the  results  of  similar  work 
performed  by  Ghajar  et  al.  (2008).  Their  experiment  was  selected  for  comparison  because  they 
tested  a  stainless-steel  tube  with  a  667  pm  diameter  which  is  similar  to  that  of  the  current  work. 
The  stainless-steel  tube  used  in  their  experiments  had  a  surface  roughness  of  ~2.9pm.  With 
similar  diameters,  differences  in  the  resulting /would  likely  be  due  to  differences  in  the  surface 
roughness.  The  channels  of  Ghajar  et  al.  (2008)  have  greater  roughness  (two  to  ten  times  greater) 
for  their  stainless-steel  tube  than  the  current  experiment  with  the  niekel  tube.  The  lower  surface 
roughness  would  likely  delay  the  transition  from  laminar  to  turbulent  flow.  Figure  53  shows/ 
(corrected  for  accessory  fittings)  for  the  present  experiments  and  the  work  of  Ghajar  et  al. 

(2008).  In  addition,  Figure  53  shows  the/curves  for  Hagen-Poiseuille  flow  (White,  1974). 
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Figure  53.  Comparison  of  Experimental  Friction  Factor  without  Head  Loss  to  that  of 
Ghajar  et.  al  (2008),  Hagen-Poiseuille  Flow  (White,  1974) ,  and  Blausius  (White,  1974). 


Many  experiments  show  departure  (yellow  eurve)  near  a  Reo  of  1000  whieh  continues  to 
increase  until  transition  to  turbulent  flow.  Few  tube  experiments  show  departure  from  the  Hagen- 
Poiseuille  curve  so  close  to  transition  as  does  the  Ghajar  et  al.  (2008)  data. 


The  experimental  data  curve  for  our  circular  channels  appears  to  depart  from  Darcy  friction  near 
1 100  Rcd  (Figure  52).  For  validation,  correlations  provided  by  Samoilenko  (1966)  were  used  to 
estimate  the  flow  departure  from  the  Hagen-Poiseuille  curve  and  the  beginning  and  end  of  the 
transition  from  laminar  to  turbulent  flow.  The  equation  for  departure  from  Hagen-Poiseuille  flow 
(smooth  tube)  is 


/  0.0065  \ 

Rbq  —  754  exp  I  — ^ —  1  (21) 

with  s/D  greater  than  0.007  (Idelchik,  1986).  Because  the  surface  roughness  is  unknown  for  the 
present  channels,  the  equation  is  solved  for  s  using  Reo  estimated  from  the  point  of  departure 
from  the  Hagen-Poiseuille  curve.  Figure  54  shows  the  experimental  data  (without  the  head  loss 
due  to  plumbing  accessories)  with  Hagen-Poiseuille  flow  (Darcy  friction  factor),  the  Blasius 
curve  for  turbulent  flow  (White,  1974),  and  a  fit  from  Ghajar  et  al.  (2008).  The  Samoilenko 
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correlation  for  departure  (Equation  3)  results  in  a  roughness  of  2.4|um  when  using  a  Reo  of 
-1000.  If  8  is  2.4|um,  then  s/D  is  less  than  0.007,  and  the  correlations  are  not  reliable  there.  The 
roughness  estimation  of  2.4  qm  is  far  from  the  expected  value.  If  the  channels  are  as  smooth  as 
expected,  such  early  departure  from  laminar  flow  is  not  expected.  This  would  mean  that  the  point 
at  which  the  experimental  data  crosses  the  Hagen-Poiseuille  curve  (White,  1974)  is  mainly  due 
to  experimental  uncertainties.  Using  a  log  plot  of  the  experimental  data  (Figure  54),  there  is  not  a 
signifieant  profile  change  near  the  crossing  of  the  Hagen-Poiseuille  flow  curve.  The  only 
significant  profile  change  occurs  near  3600  Reynolds.  This  means  little  departure  before 
transition  and  confirms  the  data  of  Ghajar  et  al.  (2008).  Thus,  the/profile  of  the  current 
experiments  should  more  closely  follow  the  Hagen-Poiseuille  curve. 


Figure  54.  Significant  Change  in  the  Flow  Near  a  Reo  of  3600  Indicated  by  the  Upward 

Bend  in  the  Friction  Factor  Curve 


Using  the  Samoilenko  correlations  again  but  for  the  start  of  transition  from  laminar  to  turbulent 
flow  (Idelchik,  1986), 


Re^  =  1160 


^e/D^ 


0.11 


(22) 


61 


The  new  estimated  roughness  is  a  very  small,  20  nm,  which  is  close  to  the  estimated  15  nm  of 
Papautsky  (1999).  When  Equation  4  is  used  with  the  Ghajar  et  al.  (2008)  data,  Rci  is  -2100  for 
the  given  surface  roughness  of  2.9  |um.  For  the  current  experimental  channels,  the  final  transition 
to  fully  turbulent  flow  is  estimated  as  (not  recorded  in  present  experimental  data) 


Re2  =  2090 


f  aO.0635 


(23) 


The  Re2  value  is  4030  with  a  corresponding /value  of  0.0348  (Idelchik,  1986).  This  agrees  well 
with  the  Moody  diagram  for  smooth  channels.  The  experimental  data  profile  slope  is  lower 
(Figure  54)  compared  to  Hagen-Poiseuille  flow  (Darcy  friction)  but  is  not  considered  critical  for 
this  analysis.  The  short  length  of  the  channels,  sharp  channels  inlets,  and  smooth  channels  could 
affect  the  flow  characteristics  and,  thus,  cause  a  difference  between  the  theoretical  and 
experimental  values.  This  section  has  shown  that  the  channels  have  a  low  pressure  drop  which  is 
important  for  aircraft  applications. 
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7.0  CONTACT  RESISTANCE 


2 

One  goal  of  the  present  work  was  to  pass  100  W/cm  through  the  8mm  x  8mm  eopper  stem.  This 
task  was  to  be  done  without  the  interface  temperature  (stem  to  heat  exchanger  contact  surface) 
reaching  200‘’C.  The  experimental  stem-heat  exchanger  results  showed  that  this  was  an 
achievable  task  without  any  resistance  (any  material  samples)  at  the  interface.  Acceptable 
contact  resistance  would  keep  the  interface  temperature  low  while  heat  is  transferred  across  the 
contact  point  between  the  stem  and  the  heat  exchanger. 


In  the  current  work,  a  silver-silicon  thermal  grease  (Arctic  Silver  5)  was  used.  Results  in  the 
Convective  Heat  Transfer  Section  reveal  that  contact  resistance  had  some  influence  on  the 
experiments.  The  contact  resistance  issues  were  due  to  daily  mechanical  positioning  and  setup 
(jumps  in  plotted  data).  Detailed  work  was  not  planned  or  conducted  for  finding  the  contact 
resistance  for  this  initial  phase  of  the  overall  project. 


Regardless,  a  relatively  simple  analysis  of  the  contact  resistance  was  performed  to  aid  in  future 
work.  Testing  was  performed  on  the  stand  with  and  without  grease  in  the  system  to  characterize 
the  thermal  behavior  of  the  thermal  grease.  Figure  55  shows  a  plot  of  the  calculated  heat  transfer 
to  the  water,  and  Figure  56  shows  the  calculated  heat  flux  through  the  stem.  The  experimental 
conditions  that  are  assessed  are  for  a  mass  flow  rate  of  7g/s  flow. 


Flow  Rate  g/s 

Figure  55.  Rate  of  Enthalpy  Rise  for  the  Cooling  Fluid  with  and  without  Thermal  Grease 
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Conduction  Rate  through  Copper  Stem  (W 


Figure  56.  Heat  Transfer  Rate  through  the  Copper  Stem  with  and  without  Thermal  Grease 
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8.0  CONCLUSIONS 


A  heat  exchanger  with  mini-channels  was  designed,  tested  and  met  a  requirement  to  transfer  100 
W/cm^  of  thermal  energy  to  the  cooling  fluid  with  a  surface  temperature  less  than  200°C.  The 
heat  exchanger  performed  well,  and  it  is  expected  that  it  can  perform  sufficiently  with  resistance 
loads  placed  in  the  interface  region  (new  module  materials).  The  two  goals  of  100  W/cm  and 
200°C  were  met  with  the  testing  performed,  and  a  sufficient  estimation  of  the  pressure  loss  was 
obtained  for  the  channels. 


The  test  stand  and  heat  exchanger  showed  relatively  good  effectiveness  in  transferring  thermal 
energy.  The  heat  input  from  the  heater  was  90  W,  and  the  heat  transfer  through  the  stem  into  the 
heat  exchanger  was  71-75  W.  The  flux  through  the  stem  was  110-117  W/cm^  with  a  maximum 
stem  temperature  of  1 12°C  at  the  interface.  The  Nud  values  for  the  flow  range  were  9-23  which 
followed  several  different  heat  transfer  correlations  reasonably  well.  The  flow  characterization 
showed  that  flows  up  to  a  Rcd  of  3600  agree  well  with  the  theoretical  pressure  loss.  The  pressure 
increase  varies  gradually  with  Rcd  up  to  the  initiation  of  turbulent  flow  with  only  an  18.5  kPa 
increase  over  the  entire  laminar  flow  range  (Rcd  up  to  3600). 


The  transition  to  turbulent  flow  at  a  Rcd  of  3600  for  the  present  experiments  is  greater  (by  1000) 
than  other  works  reviewed  and  is  likely  due  to  the  very  low  surface  roughness  of  the  channels 
(0.2  pm).  Although  there  is  a  high  uncertainty  in  the  friction  factor  value,  the  uncertainty  in  the 
Reynolds  value  is  not  great  enough  to  address  the  difference. 


With  correction  of  the  repeatability  issues  due  to  the  load  cell  alignment,  the  system  is  ready  for 
the  second  phase  of  the  project.  This  second  phase  involves  the  testing  and  property  evaluation 
of  new  materials  that  will  be  studied  for  potential  use  in  SiC  devices.  The  test  stand  will  provide 
an  accurate  way  to  determine  temperatures  at  the  material  interfaces. 


Eventually,  the  project  will  reach  the  third  phase  where  actual  modules  will  tested.  Testing  with 
the  designed  heat  exchanger  is  expected  to  continue.  The  length  of  the  channels  needs  to  be 
considered  before  continuing  with  real  module  testing  as  only  one  length  was  used  in  the 
experiment.  Another  consideration  that  should  be  made  is  the  inlet  and  exit  connections.  A 
probable  design  would  be  a  triangular  plenum  with  internal  channel  dividers.  The  plenum  should 
keep  the  area  constant  by  varying  the  divider  size  from  the  inlet  tube  to  the  channels  to  minimize 
pressure  loss.  Lastly,  CFD  simulations  provided  insight  into  the  effect  of  the  heated  length  on 
Nusselt  number  estimates  and  flow  and  convective  heat  transfer  behavior.  CFD  simulations 
should  also  be  performed  concurrently  with  future  experiments  to  provide  additional 
understanding. 
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APPENDIX  A 

SAMPLE  CALCULATIONS 


This  appendix  shows  how  various  quantities  were  caleulated  using  seleeted  test  eonditions. 

A I  Convective  Heat  Transfer 

To  ealculate  the  rise  in  enthalpy  the  following  equation  is  used 

^in  *  ^in  T  Q  ~  ^out  *  ^out  •  (A-^^ 

Using  steam  tables  for  the  enthalpy  values  at  an  inlet  temperature  of  20°C  and  an  exit 
temperature  of22AS°C  (Sonntag,  1991),  the  heat  transferred  to  the  eleven  mini-channels  is 

0.007 in  *  83940i„  +  Q=  0.007 ^ut  ^  *  94320^^^  ^  (A.  2) 

72.66  W  (A.  3) 


A2  Conduction  Heat  Transfer 

An  energy  balance  is  important  for  verification  of  the  energy  transferred  to  the  water.  The  stem 
contains  thermocouples  that  are  used  to  calculate  the  rate  of  thermal  energy  transfer  to  the  heat 
exchanger.  The  equation  describing  one-dimensional  axial  conduction  through  the  stem  is 

AT 

=  (A  4: 

Heat  transfer  through  the  stem  (for  7  g/s  fluid  flow  in  the  heat  exchanger)  with  a  temperature 
drop  of  13.14“C,  a  conductivity  of  the  copper,  k,  of  402  W/mK,  a  heat  transfer  area.  A,  of  0.008^ 
m  ,  and  a  differential  length.  Ax,  of  0.004749  m  gives 


Qs 


=  402  *  0.008^  * 


13.14 

0.004749 


(A.  5) 


The  conduction  heat  transfer  rate  through  the  stem  is  7 1 .2  W  and  the  enthalpy  rise  in  the  water  is 
72.66  W.  This  shows  good  agreement,  less  than  a  2%  difference. 
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A3  Heat  Loss 


Using  a  one-dimensional  heat  eonduetion  model,  the  heat  loss  through  the  heater  body  and  base 
are  evaluated  here.  Figures  A-1  and  A-2  show  the  physieal  situations  that  are  analyzed  in  this 
section. 


Copper  Heater 
Body 


Copper  Heater 
Base 


Qside 


Insulation 


0.0508  m 


T 


Q  base 


Ri  =  0.01207  m 
82  =  0.009525  m 


Figure  A-1.  Components  of  the  Copper  Block  Containing  the  Heater 


Insulation 


Anodized  Aluminum 


Figure  A-2.  Components  helow  the  Copper  Heater  Body 


Values  required  for  the  heat  loss  calculations  (for  a  heat  exchanger  mass  flow  rate  of  7  g/s)  are; 
air  temperatures  of  44.5  and  36.52‘’C  (lower  and  upper  regions  for  an  average  of  40.5“C),  and 
copper  body  temperatures  of  123.45  and  105.74°C  (lower  and  upper  regions  for  an  average 
temperature  of  1  Id.h^C).  The  outer  temperatures  of  the  insulation  were  measured  to  be  89  and 
7 1 .75°C  (lower  and  upper  regions  for  average  of  80.4°C).  The  temperature  at  the  base  is 
123.45‘’C.  The  heat  sink  for  the  heater  base  thermal  energy  loss  (air  and  table)  is  20°C. 
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In  terms  of  the  thermal  resistances  for  the  insulation  and  air  convection,  the  heat  transfer  rate  can 
be  represented  as 


Q 


ATi  _  ATa 
Rj  Ra 


(^.6) 


ATi  and  ^Ta  are  the  temperature  differences  (°C)  across  the  insulation  and  from  the  insulation 
exterior  to  the  air,  respectively.  The  thermal  resistances,  Ri  and  Ra,  are  for  the  insulation  and  the 
air.  The  heat  transfer  rate  can  also  be  expressed  in  terms  of  the  thermal  resistances  in  series 


(Tn  -  Ta) 
Ri  +  Ra 


(^.7) 


where  the  inner  insulation  temperature  is  Tn.  Here,  the  convective  heat  transfer  coefficient,  h, 
was  estimated  using  a  relation  for  free  convection  from  a  vertical  surface  (from  Incropera  and 
DeWitt  (eq.  9.26  and  pg.  493))  taken  as  7.3  W/m^K,  the  area  of  heat  transfer.  A,  is  2.138E-4  m^. 
Ax  is  0.00508  m,  and  the  thermal  conductivity  of  the  insulation  is  0.04  W/m  K.  In  addition,  Tn  - 
Ta  is  75°C.  With  substitution  of  the  necessary  values,  the  heat  loss  (Incropera  and  DeWitt,  2002) 
is  estimated  as 


75 

,  ^  0.01207^ 

^0.009525-*  I  1 

0.04  *  2  n  0.0508  ^  7.3  *  0.00385 


75 

18.51  +  10.38 


1.4  Watts. 


(71.8) 


The  heat  loss  through  the  body  of  the  copper  block  heater  assembly  was  approximately  1.5%  of 
the  axial  conduction  (90.74  W). 

To  calculate  the  heat  loss  through  the  copper  block  base  and  to  the  table,  the  following  equation 
can  be  used 


(Tb  -  Tt) 

"f  "f  ^A1  "f  ^A2  "f 


=  Q- 


(71.9) 


Here,  Tb  is  the  temperature  of  the  base,  and  Tt  is  the  table  temperature.  There  are  five  resistances 
to  consider,  one  for  each  component,  except  the  aluminum  which  has  two  resistances  because  of 
an  odd-shaped  cutout.  The  required  thermal  conductivities  are  0.04  W/mK  for  the  insulation, 

1.46  W/mK  for  the  Macor,  250  W/mK  for  the  aluminum,  and  4  W/mK  for  the  phenolic. 
Substituting  the  appropriate  values.  Equation  A.9  becomes 


103.5 


0.0109  +  6  +  0.094  +  0.0344  +  1.039 


=  Q  =  14.4  Watts. 


(A.  10) 


Subtracting  the  total  heat  losses  from  the  heater  output  leaves  a  net  heat  transfer  rate  of  74.26  W 
for  this  condition.  This  net  heat  transfer  rate  agrees  well  with  the  enthalpy  rise  in  the  fluid  (72.66 
W). 
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A4  Experiment  and  Correlation  Nusselt 

The  heat  transfer  to  a  heat  exehanger  ehannel  is  important  for  understanding  the  effeetiveness  of 
the  heat  exehanger  design.  Conveetive  heat  transfer  is  eommonly  deseribed  in  terms  of  the 
dimensionless  Nusselt  (Nu)  number.  Here,  Nud  is 


where 


(A.  11) 


QIA 

(.Ts-TmY 


(A.  12) 


2 

The  surfaee  area  of  the  eleven  ehannels  is  2.3957E-4  m  (D'n:L*#ehannels).  L  is  the  heated 
length  of  the  ehannels.  Here,  the  temperatures  Tm  (bulk  fluid  average  temperature)  and  Ts  (inner 
wall  temperature  estimated  from  the  stem  temperature)  are  21.8°C  and  45.3“C,  respeetively.  The 
net  heat  transfer  value,  Q,  is  72.66  W  (enthalpy  rise).  Substitution  of  the  numerieal  values  in 
Eqn.  A.  12  gives 


72.66 

,  ^  2.3957E  -  4 
(45.3  -  22.6) 


303293.4 

22.64 


W 

13396.35^—  . 
m‘^K 


Nud  is  then 


{A.  13) 


Nud 


hD 

k 


133396.35  *  0.000635 
0.602 


14.13  . 


(A.  14) 


In  the  eonveetion  discussion,  heat  transfer  correlations  are  compared  to  the  experimental  results. 
Eor  example,  Eqn  A.  15  is  for  uniform  heat  flux  and  thermal  entry  flow  conditions  where  location 
X  begins  at  the  tube  inlet: 


NUd—x  I 


3.302x*-i/^  -  1.00  forx*  <  0.00005 
1.302x*-i/^  -  0.5,  for  0.00005  <  x*  <  0.0015 
4.364  +  8.68(10^x*  exp(-41x*)  ,  for  x*  >  0.001  ■ 

xfD 


(A.  15) 


X’ 


( 


RenPr 


;) 


The  channels  tested  are  heated  at  a  central  location  along  the  flow  path,  and  the  heat  transfer 
length  varies  with  the  flow  due  to  heat  diffusion  in  the  heated  wall.  The  correlation  above  gives 
the  local  Nu  values  along  the  tube  length.  The  second  x  position  along  the  heat  transfer  length  is 
evaluated  as  an  example.  The  x  position  is  0.00109  m  from  the  start  of  the  heated  length.  With 
substitution  for  the  x  position,  the  correlation  value  is 
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Nud-x  =  1.302  *  0.0002061-1/^  _  q  5  ^  21.54 


(^.  16) 


A5  Pressure  Loss  Z-Shaped  Elbow 

The  calculation  of  the  head  loss  coefficients  is  described  in  Idelchick  (1986).  As  an  example,  the 
equation  for  non-uniform  inlet  flow  in  an  elbow  is 

Un  =  k*Cloc*Cavg  +  ^fr  (A.17) 

where  ^  is  the  pressure  loss  coefficient  (loc  refers  to  a  local  loss,  fr  is  a  frictional  loss,  and  avg  is 
a  loss  calculated  from  the  averaging  of  the  velocity  profile).  The  constant  k  depends  on  the 
velocity  profile  and  is  detailed  in  Idelchik  (1986).  For  the  current  situation  considered,  k  is  1.05. 

^"loc,  is  a  combination  of  additional  factors  (Idelchik,  1986). 

^'loc  ~  ^A^Re^loc  13) 

kA  is  a  geometrical  constant,  and  kRe  is  a  constant  based  on  the  Reynolds  number  and  the  elbow 
bend  radius,  diameter,  and  width.  Here,  the  value  for  the  two  constants  is  1 .  The  constant  ^"loc  is 
also  a  combination  of  three  more  constants  (Idelchik,  1986). 

9ioc  —  (A.  19) 

Ai  is  for  the  bend  angle  and  has  a  value  of  1 .  Bi  constant  is  obtained  from  the  bend  radius, 
diameter,  and  elbow  width  and  is  1.18  here.  Ci  is  obtained  from  the  elbow  width  and  height  and 
is  1 .  Qvg,  is  a  geometrical  constant  calculated  from  the  radius,  diameter,  width,  and  elbow  length. 
Here,  it  is  1.04.  Finally,  the  overall  value  for  ^"loc  is  1.288. 


A6  Friction  Factor:  Channels 

The  friction  factor  is  important  in  determining  the  flow  behavior  through  passages.  The  equation 
for/ is  (Wilcox,  2000) 


/  = 


(A.  20) 


AP  is  the  pressure  loss  in  the  channels,  L  is  the  channel  length,  D  is  the  channel  diameter,  p  is 
the  fluid  density,  and  w  is  the  local  average  velocity  of  the  flow  within  the  channels.  For  a  flow 
rate  of  7  g/s  (total  all  channels)  for  representative  conditions  used  in  the  experiments,  the 
velocity  is  2.03176  m/s  (each  channel),  the  density  is  989  kg/m^,  and  the  pressure  drop  was  4194 
Pa.  The  friction  factor  is  then 
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4194 


/=  f  .02159  \r 
V0.000635J  V 


989  *2.03176^^  =  0.06  (71.21) 


) 


A7  Critical  Radius:  Insulation 

An  analysis  was  performed  to  determine  how  mueh  insulation  if  any  would  be  needed  around  the 
heater  and  stem. 


Figure  A-3.  Cross-Sectional  View  of  the  Heater  Showing  the  Heater  Radius,  Outer  Radius 
of  the  Insulation,  and  Centerline  Radius  of  the  Insulation 

The  critical  radius  is  the  radius  below  which  adding  insulation  will  increase  the 
convective  heat  loss  because  the  resistance  to  heat  transfer  decreases  more  due  to  the 
increase  in  area  than  it  increases  due  to  the  insulation.  For  this  cylindrical  shape,  it  is 


r*  —  —  —  critical  radius 
h 


(A.  22) 


The  critical  radius  is  a  function  of  the  insulation  thermal  conductivity  and  the  convective  heat 
transfer  coefficient  h.  The  critical  radius  here  is  0.04/7.3  or  0.0055  m.  The  outer  radius  of  the 
heater  body  is  0.0095  m  so  adding  insulation  will  decrease  heat  transfer. 
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APPENDIX  B 

UNCERTAINTY  ANALYSIS 


BI  Parameter  Uncertainty 

Evaluation  of  the  uncertainty  in  calculations  for  the  flow  enthalpy  rise,  stem  conduction, 
Reynolds  number,  friction  factor,  and  Nusselt  number  are  described  in  this  section.  Uncertainties 
in  the  pressure,  temperature,  and  mass  flow  rate  are  described  in  later  appendices  which  also 
explain  the  calibration  methods.  Indications  of  the  uncertainties  for  length  dimensions  and 
properties  are  listed  in  Table  B-1. 

Table  B-I.  Uncertainties  for  Miscellaneous  Properties  and  Dimensions 


Length 

0.0000127 

meter 

Diameter 

0.0000127 

meter 

Density 

1%  Property  Table  Value 

kg/m^ 

Dynamic  Viscosity 

1%  Property  Table  Value 

kg/(ms) 

Thermal  Conductivity 

1%  Property  Table  Value 

W/mK 

To  evaluate  the  uncertainty,  the  following  equation  is  used  (Holman,  2001) 


Wr  = 


where  Xi  is  the  independent  variable,  Wi  is  the  individual  uncertainty,  and  R  is  the  function  being 
evaluated.  Inspection  of  Table  B-1  shows  that  some  of  the  uncertainties  are  variable  with  flow. 
Thus,  there  is  not  a  constant  wr  for  many  of  the  uncertainties.  One  evaluation  of  uncertainty  is 
shown  for  each  the  parameters  listed  at  the  beginning  of  this  section. 


B2  Stem  Conduction 

The  first  parameter  to  evaluate  is  the  one-dimensional  conduction  heat  transfer  through  the  stem 
which  is  given  by: 


AT 

Q  =  kA-  (8.2) 

k.  A,  AT,  and  Ax  which  all  have  uncertainties  which  need  consideration.  Performing  an 
uncertainty  analysis  gives 
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+ 


(5.3) 


Wq  = 


dQ 


dQ 


dQ 


\dxj^  J  \dxji 
With  substitutions,  the  above  equation  beeomes 


1/2 


Wq  = 


fAAT  (kAT  W  (kA  \ 


"  ikAAT 
+  I  ■  ..  w, 


V  Ax'^ 


Ax 


(5.4) 


The  evaluation  is  for  the  7  g/s  mass  flow  rate  eondition.  The  variables  are  the  temperature  drop 
between  two  thermoeouples  on  the  stem  (dT=13.14“C),  the  distanee  between  two  thermoeouples 
(dx  =  0.004749  m),  the  thermal  conductivity  (k  =  402  W/mK)  for  the  copper  stem,  and  the  area 
of  stem  (0.008^  m^). 


Wr 


0.008^  *  13.14 
0.004749 


(.01*402)  + 


^402  *  13.14 
,  0.004749 


(3.25  xlO-^^) 


+ 


^402  *  0.008^ 
,  0.004749 


■0.17  + 


402  *  0.008^  *  13.14 
0.0047492 


(1.27x10“^) 


(5.5) 


The  uncertainty  for  the  7  g/s  (total  flow  all  channels)  condition  is 

Wq  =  (0.527  +  1.31e  -  05  +  0.848  +  0.036)1/^  ==  +  1.2  W  .  (5.5) 


B3  Enthalpy  Rise  in  Cooling  Fluid 

The  equation  for  the  enthalpy  rise  of  the  water  in  the  heat  exchanger  is 


Qf  =  rhCpAT 


(5.7) 


where  Cp  is  the  constant  pressure  specific  heat  and  has  an  individual  uncertainty  of  0.5%.  The 
uncertainty  is  0.1 6“C  for  the  fluid  thermocouples.  The  mass  flow  rate,  rh,  has  an  uncertainty  of 
0.071  g/s.  The  overall  uncertainty  for  Qf  is 


^Qf 

Qf 


(5.8) 


where  ai  is  the  power  of  the  variable,  Wxi  is  the  individual  uncertainty,  Xi  is  the  variable,  and  the 
denominator,  Qf,  is  variable  with  flow  rate.  Expansion  of  Eqn  B.6  yields 
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^Qf 

Qf 


^cp^cpy 

Xfii  ^cp  ^T 


Upon  substitution  of  numerical  values 


(5.9) 


^Qf 

72.66 


1  *  0.0715  -  3  ,  1  *  4179  *  0.005  ,  1  *  0.16 

( - 555^ - +  ( - 77^5 - +  (^77^)^  (5- 10) 


0.007 


4179 


2.46 


The  uncertainty  changes  with  the  flow  rate  which  also  influences  the  temperature  increase  in  the 
fluid.  The  value  for  dT  at  a  mass  flow  of  7  g/s  (total  flow  all  channels)  is  2.46°C  and  the  Qf  value 
is  72.66  Watts.  The  resulting  fractional  uncertainty  is 

Wnf 

— ^  =  [0.000103  +  0.000025  +  4.235  -  3]^/^  =  6.7%  (5. 11) 

72.66 

Thus,  the  uncertainty  for  the  heat  transfer  is  approximately 


WQf  =  +4.8  W. 


(5.12) 


B4  Friction  Factor 

The  friction  factor  relation  is 


(5. 13) 


Expressing  the  velocity  (V)  in  terms  of  the  mass  flow  rate  gives 


2A55 


/  = 


Lp 


n2 


m 


(  £,2^) 


TT^  APD^p 
8  Lrh^ 


(5. 14) 


for  which  AP,  diameter,  density,  length,  and  mass  flow  have  uncertainties.  With  substitution  in 
the  uncertainty  equation 
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TT^ 


Wf  = 


yLrh^ 


+ 


'2*  APD^p 
Lm? 


Wap  )  + 


nl/2 


'S*APD^p 
,  Lm? 


Wd  + 


/APD'^ 
I  Lrh^ 


■w„ 


+ 


fAPD^p 
i  L?-rh?- 


-Wl 


■Wr, 


{B.  15) 


For  the  flow  conditions  associated  with  7  g/s  (or  0.000636  g/s  for  a  single  channel),  the  variable 
values  are  4184  Pa  for  AP,  0.000635  m  for  D,  989  kg/m^  for  p,  and  0.02159  m  for  L.  The 
uncertainty  values  are  345  Pa  for  pressure,  1.27*10'^  m  for  D  and  L,  0.0 Ip  for  p,  and  0.000071 
kg/s  for  the  mass  flow  rate.  The  result  for  this  evaluation  is  0.01477  or  a  25%  uncertainty  (for  a 
nominal  f  of  0.06). 


B5  Reynolds  Number 

The  Reynolds  number  written  in  terms  of  the  mass  flow  rate  is 

4rh 

Re  =  -  (B.  16) 

nDp  ^  ^ 

where  mass  flow  rate,  diameter,  and  viscosity  have  uncertainties.  The  individual  uncertainties  are 
0.000071  kg/s  for  the  mass  flow  rate,  1.27  x  10'^  m  for  the  diameter,  and  1%  for  the  viscosity. 
The  uncertainty  equation  is 


WRe  = 

Evaluated  partial  derivatives  give 

4 

WRe  =  - 


dRe 

dXjji 


w. 


dRe  ^ 


1/2 


(5. 17) 


TT 


2  •  21^/2 
m  \  /  m  '  ' 


(5. 18) 


The  resulting  uncertainty  for  Rcd  is 

/  1  0.007  ( 

- t(0.000071  * - )  +  - 

V0.000635  *  1020x10-6 H  J J  Iq 

r0.007> 


WRe  = 


0.007 

(  II  ) 


■1.27x10-6 


+ 


/'O.OO/^ 

I  11  7 


0.000635  *  (1020x10-6)2 


.0006352  ,,  1020x10-6 

2-,1/2 

(0.01  *  1020x10-6)^  (5.19) 


WRe  =  22 


(5.20) 
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The  uncertainty  in  Re  at  a  flow  rate  of  7  g/s  (total  flow  all  channels)  is  roughly  +  30. 


B6  Nusselt  Number 

The  convective  heat  transfer  is  commonly  described  in  terms  of  the  dimensionless  Nusselt  (Nu) 
number. 


hD 

Nud  =  — 

(B.21) 

with 

h 

(T,-T^y 

(B.22) 

The  quantity  Q  is  the  total  heat  transferred  to  the  channels  (expected  '-80W),  A  is  the  channel 
surface  area  (11  channels),  Tg  is  the  channel  surface  temperature,  and  Tm  is  the  mean  bulk 
temperature.  Combining  Equations  B.14  and  B.15 

Q/A  ^ 

Nud  =  ^ 

(B.23) 

Q  *  D 

A*k*(Ts-  Trn) 

(S. 24) 

Q  *  D 

~  DnL  *ll*k*(Ts-  Trn) 

(B.25) 

-u  -  ^ 

nL*ll*k*(T,-Tm) 

(B.26) 

The  uncertainty  in  Q  (physical)  is  6  W,  in  L  it  is  1.27*10'^m,  in  k  it  is 
(Ts-Tm)  it  is  0.16  °C.  Using  Eqn.  B.l 

1%  of  0.56  W/mK,  and  in 

1  T/  1  \^  (  Q 

(  Q  f  Q 

Using  the  7  g/s  heat  transfer  length  for  E,  1 1  mm. 
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Wnu  = 

+ 

*  0 


lln 


0.011  *0.56  *22.64 
80 


80 


0112  *  0.56  *  22.64 


1.27  *  10"^ 


0.562  *0.011*22.64 
2  '  80 


0.01 


■“)  +(o: 


011  *  0.56  *22.642 


2-|1/2 

).16j 


(5.28) 


Wwu~l-3 


(5.29) 
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APPENDIX  C 

CALIBRATION  METHODS 


The  thermocouples,  pressure  transducers,  and  the  flow  meter  were  calibrated.  A  load  cell  is 
present  in  the  system  and  typically  would  be  calibrated.  However  since  the  load  cell  had  a 
calibration  certiflcate,  it  was  not  calibrated.  The  platinum  resistance  thermometer  (PRT)  (used  as 
a  calibrating  unit  for  the  thermocouples)  calibration  was  performed  by  an  in  house  calibration 
department.  In  addition,  two  points  were  verified  in  the  RZPS  Thermal  Laboratory:  the  freeze 
point  and  the  normal  boiling  point  for  a  distilled  water  sample.  All  curve  fit  equations  were 
obtained  using  KaleidaGraph  software. 


Cl  Thermocouple  Calibration 

The  equipment  for  the  thermocouple  calibration  included  a  Hart  Scientific  1502A  PRT  reader. 
Hart  Scientific  7320  cooler/bath,  and  Hart  Scientific  5628  PRT.  The  calibration  bath  contained 
Dow  Corning  200.50  silicone  oil.  The  equipment  is  shown  in  Figure  C-1.  The  indicated  bath 
temperature  was  used  as  a  reference  set  point  only.  All  calibrations  were  performed  with  the 
Agilent  unit  set  on  internal  reference  compensation.  The  Agilent  unit  was  operated  for  at  least  30 
minutes  before  any  calibrations  or  testing  was  performed. 


Recording  of  the  thermocouple  readings  started  when  the  bath  temperature  reached  the  expected 
value.  This  was  monitored  by  the  PRT,  and  the  bath  temperature  was  considered  to  be  at  steady- 
state  conditions  when  a  ±  0.1  °C  fluctuation  per  hour  was  obtained. 


The  calibration  bath  had  a  temperature  limit  of  150‘’C,  but  test  temperatures  for  the  stem  could 
possibly  reach  200°C.  Thus,  a  different  calibration  bath  was  used  for  the  higher  temperatures, 
and  it  was  capable  operating  safely  at  250°C.  Trials  with  the  bath  showed  that  linear  fits  were 
reliable  up  to  210°C. 


The  setup  was  such  that  the  thermocouples  and  the  PRT  were  approximately  at  the  center  depth 
of  the  bath  oil.  The  PRT  and  loose  thermocouples  were  attached  to  an  external  adjustable  stand. 
To  hold  the  copper  heater  off  the  bath  base,  an  insulating  block  of  Macor  was  used.  This  limited 
the  effects  of  the  steel  bath  housing  from  creating  a  conduction  path.  This  method,  however,  left 
the  PRT  body  improperly  submerged  in  the  oil.  To  compensate  for  the  improper  level,  the  entire 
PRT,  bath  oil  tank,  and  the  thermocouples  were  covered  by  a  layer  of  insulation  so  that  the 
internal  chamber  air  and  oil  would  then  be  at  the  same  temperature.  Figure  C-1  shows  the 
insulation  that  covered  the  bath.  The  thermocouples  were  calibrated  in  two  groups,  and  the  same 
equipment  and  method  were  used  for  each  group. 


The  thermocouple  calibrations  were  performed  with  equipment  that  would  be  present  during  the 
actual  experiments  (i.e.,  wire  extensions,  data  acquisition  equipment,  and  connections).  There  are 
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six  thermocouples  located  in  the  heater  stem  made  of  copper.  These  thermocouples  were 
calibrated  after  they  were  mounted  to  the  stem. 


Figure  C-1.  Photo  of  Insulation  Used  to  Cover  the  Bath  and  Calibration  Components  to 
Reduce  Errors  Caused  hy  Conduction  through  the  PRT  Body 


Cl.l  Fluid  Thermocouples 

The  fluid  thermocouples  (TC  #  13  and  14)  at  the  inlet/exit  of  the  heat  exchanger  have  (1/16  in 
diameter)  stainless-steel  sheaths.  To  allow  for  quick  change  and  to  lower  sheath  conduction 
effects,  open  tipped  probes  were  chosen.  Calibration  of  the  fluid  thermocouples  required  a  small 
AT  between  the  two  thermocouples  for  a  given  PRT  reading.  Since  the  enthalpy  rise  was  used  as 
the  energy  balance  check,  the  temperature  increase  (AT)  was  critical.  The  fluid  thermocouples 
were  calibrated  for  a  small  temperature  range  since  the  expected  fluid  temperatures  were 
expected  to  be  in  the  range  30-35°C.  At  increments  of  3  °C,  the  thermocouples  were  calibrated 
from  14  to  39°C. 


Various  curves  were  fit  to  the  thermocouple  data.  The  curves  are  listed  in  Table  C-1.  The 
residuals  of  the  fit  curves  were  analyzed,  and  the  best  fit  determined.  The  best  fit  was  determined 
to  be  a  linear  curve  with  the  standard  deviation  between  them  being  0.04‘’C.  The  goodness  of  fit, 
R  ,  is  0.99999  and  1  for  each  fit.  The  equations  for  the  two  thermocouples  are  (TC  13  inlet  fluid) 
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T,n  =  0.99199  *  X  +  0.59302 


(C.l) 


and  (TC  14  exit  fluid) 


TexU  =  0.98815  *  X  +  0.73129. 


(C.2) 


The  X  value  is  the  voltage  reading  from  the  Agilent  unit.  The  Agilent  unit  is  set  for  T-type 
thermocouples  which  automatically  fits  the  voltage  readings  to  the  standard  thermocouple  curve 
with  the  internal  electronic  ice  point. 


Cl. 2  Stem  and  Miscellaneous  Thermocouples 

The  stem  thermocouples  are  0.005  in  butt  welded  T-type.  The  thermocouples  are  glued  into 
0.020  in  diameter,  0.030  in  deep  holes  (Figure  C-2).  The  thermocouples  were  secured  with 
Omegabond  1500  °F  ceramic  epoxy  before  calibration.  Calibrating  the  thermocouples  in  place 
eliminated  re-calibration  if  the  thermocouples  were  broken  during  installation.  Calibration 
settings  were  “random”  values  between  20°C  and  142°C  for  the  all  thermocouples  that  would 
experience  temperatures  above  40“C  (stem,  chamber,  and  body). 


Figure  C-2.  Stem  with  Thermocouple  Holes 

Various  calibration  curves  were  also  fit  for  these  thermocouples,  and  the  best  fit  was  linear 
(Table  C-1).  Except  for  first  and  last  setting,  the  standard  deviation  from  TCI  to  the  other  five 
TC"s(all  stem  thermocouples)  is  approximately  0.05“C.  The  differences  among  the  readings  for 
thermocouples  1-6  were  most  important  (TC  7-12  are  for  external  readings).  The  lowest  value 
was  0.99999. 
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Table  C-1.  Curve  Fits  for  Thermocouples  1-12  (Stem,  Body,  and  Chamber) 


TC#  Linear  Fit 

1  y  =  0.9992X  +  0.0748 

2  y  =  0.9991x  +  0.0489 

3  y  =  0.9986x  +  0.0519 

4  y  =  0.9992X  +  0.002 

5  y=  1.0066X-  0.2135 

6  y=  1.007X-  0.2154 

7  y  =  1 .0057X  -  0.046 

8  y  =  1 .0059X  -  0.0046 

9  y  =  1 .0064X  +  0.0087 

10  y=  1.0049X  +  0.1155 

11  y  =  0.9982X  +  0.2984 

12  y  =  0.9976X  +  0.3252 


The  temperature  uneertainty  due  to  the  voltage  resolution  of  the  Agilent  unit  is  O.loC.  Using  the 
twiee  the  standard  deviation  of  a  nominal  measurement  during  ealibration,  the  uneertainty  due  to 
ealibration  error,  weal,  is  estimated  to  be  0.08  oC  for  the  fluid  thermoeouples  and  O.loC  for  the 
stem  thermoeouples.  The  random  error  due  to  voltage  fluetuation  is  estimated  to  be  O.loC.  To 
estimate  the  uneertainty  in  the  thermoeouple  measurements,  the  contributions  to  uncertainty  from 
random  error,  wrand,  instrumentation  error,  winst,  and  calibration  error,  weal,  are  combined: 

W  —  (Wrand  Winst  Weal  )  (C-  3) 


Thus,  the  overall  uncertainty  in  the  temperature  for  the  fluid  and  stem  thermocouples  is  0.1 6“C 
and  0.17°C,  respectively. 


C2  Pressure  Transducer  Calibration 

Two  pressure  transducers  were  calibrated  (designated  INLET  and  EXIT)  using  a  portable 
pressure  calibrator  (Druck,  DPI  610).  The  unit  has  an  accuracy  of  0.025%  of  full  scale  (fs)  with 
macro  and  micro  pressure  adjustment  capabilities.  The  Druck  pressure  range  was  14.7-300  psia. 
The  transducers  (Omega  PX303-050A5V)  had  a  range  of  0-50  psig  with  an  output  voltage  of 
0.5-5  VDC. 


The  inlet  and  exit  transducers  were  calibrated  separately  using  the  Druck  unit  with  each  having 
24  random  levels  for  calibration.  Each  calibration  point  was  allowed  to  settle  for  30  seconds 
before  recording  the  voltage  output  values.  The  recorded  values  for  each  transducer  was  then 
averaged  for  each  setting  and  used  for  curve  fitting.  Eigure  C-3  is  a  schematic  of  the  setup  used 
for  calibration.  The  calibration  fits  for  the  inlet  and  outlet  pressure  transducers  are 


Pin  =  10.083  *  volts  —  19.547 


(C.4) 
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Pout  =  10.069  *  volts  -  19.568.  (C.  5) 


Druck 


Transducer  to  Calibrate 


Figure  C-3.  Schematic  of  Setup  for  Pressure  Transducer  Calibration 

The  uneertainty  in  the  pressure  measurement  due  to  the  voltage  resolution  of  the  Agilent  unit  is 
negligible  relative  to  the  other  uneertainties.  Using  the  twiee  the  standard  deviation  of  a  nominal 
measurement,  the  uneertainty  due  to  random  error,  Wrand,  is  estimated  to  be  0.007  psi  for  the 
pressure  transdueers.  The  uneertainty  in  the  pressure  measurement  due  to  ealibration  is  estimated 
to  be  0.05  psi.  To  estimate  the  overall  uneertainty  in  the  pressure  measurements,  the 
eontributions  to  uneertainty  from  the  random  error,  wrand  and  ealibration,  weal,  are  eombined; 

W  =  (Wrand^  +  Wcal^)^^^  (C.  6) 

Thus,  the  overall  uneertainty  in  the  pressure  is  estimated  to  be  +  0.05  psi  at  the  49.0  psi  level. 

C3  Flow  Meter  Calibration 

The  flow  meter  ealibrated  for  the  system  was  an  Omega  (FTB-9504)  preeision  low  flow  turbine 
meter  with  transmitter  (FLSC-61).  Before  flow  ealibration  eould  be  performed,  the  transmitter 
eleetronies  required  zeroing.  The  proeess  involved  a  frequeney  generator,  DVOM,  and  an 
oseilloseope.  The  adjustment  of  the  transmitter  was  without  issue.  Details  of  the  proeess  to  set 
the  eleetronie  zero  are  in  the  instruetion  manual  for  the  transmitter. 
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The  calibration  of  the  flow  meter  was  performed  with  the  unit  installed  in  the  test  stand  so 
potential  effects  from  other  components  are  taken  into  account.  Figure  C-4  shows  the  calibration 
setup. 


□ 

Pump 


□ 

Filter 


^  Flowmeter 


Flow 


Graduated  Cylinder 


Scale 


Figure  C-4.  Flow  Meter  Calibration  Setup  Schematic 

The  flow  in  the  system  begins  with  a  pump  whose  suction  hose  is  submerged  in  a  bucket  of 
distilled  water.  The  water  flows  from  the  pump  through  the  filter  and  heat  exchanger  to  the 
flowmeter  and  is  collected  in  a  graduated  cylinder  (1000  mL).  The  flow  was  provided  by  a 
constant  pressure  pump  (Aquatec  DDP  5800).  The  flowmeter  signal  was  monitored  by  the 
Agilent  data  acquisition  unit.  On  and  off  flow  was  controlled  by  a  quick  turn  ball  valve  located 
above  the  graduated  cylinder. 


The  flow  meter  was  calibrated  using  mass  accumulated  with  time.  Since  the  temperature  of  the 
water  entering  the  test  section  was  controlled  to  always  be  the  same  value  (and  same  density),  it 
it  reasonable  to  calibrate  the  flow  meter  to  yield  a  mass  flow  rate  rather  than  a  volumetric  flow 
rate.  Each  test  time  was  random  in  the  range  of  36  to  109  seconds  for  slow  flow  rates.  For  high 
flow  rates  (1.4-22  g/s),  the  cylinder  was  nearly  full.  The  flow  settings  used  are  random  by 
estimating  the  next  setting  as  falling  between  the  previous  flow  value  and  one  full  turn  of  the 
valve.  An  electronic  scale  was  used  to  monitor  the  fluid  mass.  The  flow  time  was  found  by 
reviewing  the  data  files  provided  by  the  LabView  program  that  reads  and  records  the  Agilent 
unit  values.  The  mass  indicated  by  the  scale  and  divided  by  the  mass  accumulation  time  gave  the 
mass  flow  rate  corresponding  to  the  voltage  reading  of  the  flow  meter. 


It  was  found  during  data  analysis  that  the  1  -2  second  period  required  for  starting  and  stopping  the 
flow  with  the  valves  provided  an  error  of  a  few  percent.  To  adjust  for  the  varied  flow  at  the  start 
and  end,  the  time  from  zero  to  full  flow  was  used  to  average  the  mass  flow  for  that  same  period. 
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The  averaged  mass  was  removed  from  the  total  flow  amount  and  the  time  was  taken  from  the 
total  time.  Smoother  eurve  fits  were  then  obtained  after  the  eorrection. 


The  ealibration  data  was  fit  with  various  order  eurves  to  find  the  best  fit,  and  a  fifth  order 
polynomial  provided  the  best  fit  for  the  mass  flow  rate: 


rh  =  0.42258  +  4.8607  *  x  —  0.16123  *  +  0.080099  *  x^ 

-0.009851  0.000444  (C.  7) 


After  the  eurve  aceuraeies  were  eheeked  with  the  manufaeturer  speeifieations  (0.1 1%  of  the  flow 
rate),  physieal  flow  eheeks  were  performed.  Prior  problems  with  flow  meter  failure  warranted 
the  cheek.  The  calculated  error  for  the  measured  flow  and  the  calibration  points  are  shown  in 
Figure  C-5.  The  relatively  larger  error  at  low  flows  is  expected  to  be  due  to  the  range  of  flow 
frequencies  chosen  for  the  transmitter.  The  selected  flow  calibration  range  was  slightly  too  great 
for  these  low  flow  rates.  The  frequency  range  chosen  was  350-1950  Hz.  The  error  was  calculated 
from  an  average  of  3-4  repeated  runs  at  2-8  g/s. 


0  2  4  6  8  10 

Flow  g/s 

Figure  C-5.  Error  Calibration  Curves  for  Voltage  to  Flow  Rate  Conversion 


The  uncertainty  due  to  the  voltage  resolution  of  the  Agilent  unit  is  negligible  relative  to  the  other 
uncertainties.  Using  the  twice  the  standard  deviation  of  a  nominal  measurement  (7  g/s),  the 
uncertainty  due  to  random  error,  Wrand,  is  estimated  to  be  0.01  g/s.  The  uncertainty  contribution 


87 


due  to  calibration  is  estimated  to  be  0.07  g/s.  To  estimate  the  uncertainty  in  the  mass  flow  rate 
measurements,  the  contributions  to  uncertainty  from  random  error,  Wrand,  and  calibration  error. 
Weal,  are  combined: 


W  =  (Wrand^  +  Wcal^)^^^ 


(C.8) 


Thus,  the  overall  uncertainty  in  the  mass  flow  rate  (at  7  g/s)  is  estimated  to  be  0.071  g/s  ('-'1.0%). 

One  signifleant  mistake  occurred  during  initial  calibration  attempts.  The  type  of  flow  meter  that 
was  used  relied  on  pulse  generation  which  was  provided  by  the  passing  of  a  metallic  object 
through  a  magnetic  fleld.  The  flrst  flow  meter  calibrations  were  performed  with  tap  water  to 
provide  an  uninterrupted  flow  of  water.  This  proved  to  be  a  bad  decision  because  the  level  of 
iron  in  the  water  was  high  enough  to  cause  failure  of  the  flow  meter.  However,  the  use  of 
distilled  water  resolved  the  issue. 
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APPENDIX  D 
SYSTEM  DIAGRAMS 


Computer  with 
LabView 


8  Power  Supply 
Heater 

9  Power  Supply 
for  Load  Cell, 
Pressure 
Transducers, 
and  Flow  meter 


5  Load  cell 

6  Shunt  -  Voltage  and  Current 

Sampling 


1  Pump/Cooler  -  Sinusoidal 
Temp 

2  Filter  -40micron 

3  Heat  Exchanger 

4  Heater  w/Cu  heat  spreader 


7  Flow  Meter 


Figure  D-I.  Block  Diagram  of  Electronic  Cooling  Components 
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APPENDIX  E 
WIRING  SCHEMATIC 


Capacitors 


'Current 


Flow  meter 


Pressure 

Transducers 


Ftt^stor 


Pot*ntion>«t«r 


Computer 


Heater 


Heater  Power 
Supply 


Shunt 


Load 

Cell 


Load  Apparatus 


Agilent  DAQ 


Thermocouples 


Power  Supply 


Heater  and  Stem 


Pressure 

Transducer 


(,0 

Voltage  Divider 


Figure  E-I.  Wiring  Map  of  Electronic  Cooling  Test  Cell 
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APPENDIX  F 

LAB  VIEW  DATA  ACQUISITION  PROGRAM 


The  LabView  program  is  located  on  computer  01408189  for  the  electronic  cooling  experiment 
and  is  backed  up  on  a  CD.  The  file  name  is  “Experiment  Runs  Agilent  34970A  Advanced 
Datalogger”.  The  LabView  program  was  developed  with  the  help  of  Victor  Tsao.  The  program 
was  developed  to  only  record.  No  component  control  is  capable  with  the  designed  program.  The 
front  screen  allows  operator  input  for  file  names  which  has  a  format  based  on  testing  parameters 
The  program  also  provides  real  time  numerical  and  graphical  outputs  viewable  by  the  operator. 
Values  such  as  power  input  to  the  heater,  thermocouple  readings,  pressures,  mass  flow  rate,  and 
load  are  available  on  the  front  screen.  The  program  also  allows  the  operator  to  adjust  the  time 
increment  for  recording  data  to  a  file.  Table  F-1  gives  the  channel  designation  numbers  for  the 
experiment  instrumentation.  Figure  F-1  is  the  front  screen  view  of  the  LabView  program. 


Table  F-I.  Channel  Locations  in  Agilent  DAQ  Unit  for  Each  Experimental 

Instrument 


101 

Heater  Stem  TC  Position  1 

111 

Heater  Body  TC  Inner  &  Upper  Position 

102 

Heater  Stem  TC  Position  2 

112 

Upper  Chamber  Air  Temperature 

Heater  Exchanger  Inlet  Fluid 

103 

Heater  Stem  TC  Position  3 

113 

Temperature 

104 

Heater  Stem  TC  Position  4 

114 

Heater  Exchanger  Outlet  Fluid  Temp 

105 

Heater  Stem  TC  Position  5 

115 

Fluid  Inlet  Pressure  Transducer 

106 

Heater  Stem  TC  Position  6 

116 

Fluid  Outlet  Pressure  Transducer 

107 

Lower  Chamber  Air  Temperature 
Heater  Body  TC  Outer  &  Lower 

117 

Cooling  Fluid  Flowmeter  Transducer 

108 

Position 

Heater  Body  TC  Outer  &  Upper 

118 

Heater  Amperage  (Shunt) 

109 

Position 

Heater  Body  TC  Inner  &  Lower 

119 

Load  Cell 

110 

Position 

120 

Heater  Voltage 
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Figure  F-1.  LabView  Front  Screen  View  of  Control  Program 
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APPENDIX  G 

COPPER  MATERIAL  QUANTIFICATION  (EDAX) 


HtAT  focoi>*ii>  ^Lt'^vcA’  Tsho  2/  G /2-®*S 


c ; \edax32\genasis\genspc . spc 

L«b«l : 

kV:25.0  Tilt: 0.0  Take-off : 35 . 0  Det  TYpe:SDTH+  Res: 152  Anp.T:6.4 

FS  :  16103  Lsec  :  30  6-Aug-2008  13:49:40 

Counts 

ML 

.2k 

Ok 
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\ 
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.50  15 

.50  20 

.50  25 
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.50  35 

. 50  keV 

EDAX  ZAF  Quantification  ( Standardlass ) 

Element  Nomallzed 
SEC  Table  :  Default 

Element  Wt  %  At  K-Ratio  Z_ A_ F 

CuK  100.00  100.00  1.0000  1.0000  1.0000  1.0000 

Total  100.00  100.00 

Element  Net  Inte.  Bkqd  Inte.  Inte.  Error_ P/B 

CuK  84.63  1.80  2.03  47.02 


Figure  G-I.  EDAX  (Energy  Dispersive  X-ray  analysis)  Quantification  of  Copper  Used  to 

Manufacture  the  Heater  and  Stem 
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NOMENCLATURE 


ID  One  -dimensional 

2D  Two-dimensional 

3D  Three-dimensional 

a  Width  of  reetangular  ehannel  (m) 

A  Area  (m^) 

AIN  Aluminum  nitride 

b  Height  of  reetangular  channel  (m) 

CFD  Computational  fluid  dynamics 

Cp  Constant  pressure  specific  heat 

(J/kgK) 

CTE  Coeff.  of  thermal  expansion  (1/°C) 

D  Diameter  (m) 

Dh  Hydraulic  diameter  (m) 

DBC  Direct  bonded  copper 

DVOM  Digital  Volt-Ohm  meter 

/  Friction  coefficient 

F  Cross-sectional  area  (m  ) 

Gz  Graetz  number  -  Re  Pr/(x/d) 

g  Gravitational  constant  (m/s  ) 

h  Convective  heat  transfer  coefficient 

(W/m^K) 

He  Micro-hardness  of  Copper  (MPa) 

HX  Heat  exchanger 

IGBT  Insulated  Gate  Bipolar  Transistor 
K  Pressure  loss  coefficient 


Kl  Pressure  loss  coefficient  -  minor 

loses 

k  Thermal  conductivity  (W/mK) 

rh  Mass  flow  rate  (kg/s) 

MOSFET  Metallic  Oxide 

Semiconductor  Field  Effect 
Transistor 

Nux  Focal  Nusselt  number 

Nud  Mean  Nusselt  number  based  on 

diameter  -  hD/k 

Nud-x  Local  Nusselt  number 

n  Area  Ratio 

P  Pressure  (Pa) 

AP  Pressure  Drop  (Pa) 

PID  Proportional  integral  derivative 

Pr  Prandtl  Number  (Cp|u/k) 

Q  Heat  transfer  rate  (W) 

Qs  Heat  transfer  rate  through  copper 

stem  (W) 

q”  Heat  flux  (W/m^) 

R  Thermal  resistance  (K/W) 

Ri  Thermal  resistance  for  insulation 

(K/W) 

Ra  Thermal  resistance  for  air  (K/W) 

Rc  Copper  thermal  resistance  (K/W) 

Rf  Thermal  resistance  fluid  channel 

Re  Reynolds  number 
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Rcd  Reynolds  number  based  on  diameter 

ipVD/vO 

Reo  Reynolds  number  for  friction  factor 
departure  from  theory 

Rci  Reynolds  number  for  beginning  of 
transition  to  turbulent  flow 

Rei  Reynolds  number  for  end  of 

transition  to  turbulent  flow 

RTD  Resistance  thermal  device 

r*  Critical  radius 

SIC  Silicon  Carbide 

Si  Silicon 

Si02  Silicon  dioxide 

T  Temperature  (“C  or  K) 

To  Inlet  temperature  (TC  1  &  4) 

Tc  Temperature  copper  (°C  or  K) 

Tin  Bulk  temperature  channel  inlet  (°C 
or  K) 

Tm  Mean  fluid  temperature  (°C  or  K) 

Tout  Bulk  temperature  channel  exit  (°C  or 

K) 

Ts  Channel  surface  temperature 

(Tw+Tm)/2  ("C  or  K) 

Ttip  Temperature  stem  tip  (°C  or  K) 

Tw  Temperature  channel  wall  (“C  or  K) 

AT  Temperature  difference  (‘’C  or  K) 

ATc  Temperature  drop  across  copper 

stem  (°C  or  K) 


ATc-f  Temperature  difference  between 
copper  and  fluid  (“C  or  K) 

TC  Thermocouple 

u  internal  energy  (J/kg) 

V  Velocity  (m/s) 

Vavg  Average  velocity  (m/s) 

Vfr  Fluid  velocity  along  friction  Surface 
(m/s) 

Vioc  Local  velocity  (m/s) 

VAC  Voltage  AC 

VDC  Voltage  DC 

X  Distance  (m) 

Ax  Change  in  distance  (m) 

Axsink  Distance  across  sink  (m) 

X*  (x/D)/Re  Pr 

a  Thermal  Epansion  Coefficient 

8  Surface  Roughness  (m) 

p  Density  (kg/m  ) 

|Uf  Dynamic  viscosity  (N/sm  ) 

Pw  Dynamic  viscosity  at  a  wall  (N/sm  ) 

^  Coefficient  of  flow  resistance 

(pressure  loss  coefficient) 
AP/(pVavg"/2) 

^"loc  Coefficient  of  local  flow  resistance 

^fr  Coefficient  of  flow  resistance  of  a 

segment  of  length  L 

Cavg  Average  coefficient  of  flow 

resistance 
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